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Abstract 

In the semiclassical approximation in which the electric charges of scalar par- 
ticles are described by Grassmann variables {Qf = 0, QiQj 7^ 0), it is possible 
to re-express the Lienard-Wiechert potentials and electric fields in the radia- 
tion gauge as phase space functions, because the difference among retarded, 
advanced, and symmetric Green functions is of order Q^. By working in the 
rest-frame instant form of dynamics, the elimination of the electromagnetic 
degrees of freedom by means of suitable second classs contraints leads to the 
identification of the Lienard-Wiechert reduced phase space containing only N 
charged particles with mutual action-at-a-distance vector and scalar poten- 
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tials. A Darboux canonical basis of the reduced phase space is found. This 
allows one to re-express the potentials for arbitrary N as a unique effective 
scalar potential containing the Coulomb potential and the complete Darwin 
one, whose 1/c^ component agrees for with the known expression. The effec- 
tive potential gives the classical analogue of all static and non-static effects 
of the one-photon exchange Feynman diagram of scalar electrodynamics. 
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I. INTRODUCTION 



Recently the rest-frame Wigner-covariant instant form of dynamics has been developed 
in Ref. [|I| for isolated systems in Minkowski spacetime M*^ starting from the case of scalar 
charged particles plus the electromagnetic field. The charges of the particles are described 
by bilinears in Grassmann variables following the scheme (called pseudo-classical mechanics) 
which uses a semiclassical approximation to quantum operators with a finite discrete spec- 
trum like spin ||^,^] , that otherwise would have no strict classical limit ; Grassmann variables 
give fermionic oscillators after quantization. The extension of this scheme to the electric 
charge is based on the experimental fact that all measurable charges are multiples of ±e, 
the electron and positron charges. Therefore, even if it is not clear whether the electric 
charge has to be considered as a quantum operator in the standard sense (except in the case 
of the existence of magnetic charges; in this case there is the Dirac quantization rule for 
the product of the electric and magnetic charges), one can consider it as a two-level system 
[which becomes a six-level system (±e, ±|e, ±|e) at the quark-lepton level] described by an 
operator with quantum e instead of h. Then one can define a semiclassical approximation 
with Grassmann variables like in the case of spin. As shown in Ref. this semiclassical 
approximation automatically implies the regularization of the Coulomb self-energies (the 
i ^ j rule). Therefore, this semiclassical approximation may be considered as an alternative 
to the extended electron models, which were introduced for regularization aims. 

The idea leading to the rest-frame instant form is to consider an arbitrary 3-1-1 splitting of 
Minkowski spacetime by means of a foliation with spacelike hypersurfaces S(r) diffeomorphic 
to . The parameter r labelling the leaves is used as a Lorentz scalar mathematical time 
parameter. For each r the leaf S(r) is defined through the embedding ^ ^(t) C M^, 
(r, a) 1-^ z^{T,a), where a are curvilinear coordinates on R^. Then one considers the La- 
grangian describing the coupling of the given isolated system to an external gravitational 
field and replaces the 4- metric with the induced metric on S(r), which is a functional of 
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z'^{t, a). In this way one gets the Lagrangian for the description of the isolated system on 
arbitrary spacehke hypersurfaces (i.e. in arbitrary accelerated reference frames in Minkowski 
spacetime) with the embedding functions z'^(r, a) as extra configuration variables describ- 
ing the hypersurface. However, there are four first class constraints at each point implying 
the independence of the description from the chosen 3+1 splitting. Thus the z'^(r, a)'s are 
gauge variables. Therefore, one can restrict the description of the isolated system to space- 
like hyperplanes T,h{t), z^{t, a) = a;^(r) + hf{T)(j^ (inertial reference frames in Minkowski 
spacetime; x^{t) is an arbitrary origin). 

Then, if one selects all the configurations of the isolated system with total timelike 4- 
momentum (they are dense in the space of all configurations), one finds that each timelike 
configuration identifies a privileged family of hyperplanes: those orthogonal to its total 4- 
momentum [Wigner hyperplanes T^wij)]- At this stage one has obtained the analogue of 
the nonrelativistic center-of-mass separation and the definition of a new instant form of 
dynamics [^], the rest-frame one [|I|. There is a decoupled point x^(r) on each Wigner 
hyperplane describing the "external" center of mass of the isolated system (thus serving as a 
decoupled point particle clock) with conjugate momentum p^. x^{t) is a canonical variable, 
but it is not covariant like the Newton- Wigner position operator: it has only covariance 
under the little group of timelike Poincare orbits. 

After the restriction to the Wigner hyperplane only four first class constraints are left: 

i) three of them say that the total 3-momentum of the isolated system vanishes (rest- 
frame condition): the natural gauge fixing for these constraints is the requirement that the 
"internal" 3-center of mass of the isolated system inside the Wigner hyperplanes coincides 
with the origin x^(r) of the coordinates a in it (see Refs. , where the group-theoretical 
results of Ref. |0] are used). In this way only "internal" relative variables describe the 
isolated system: they are either Lorentz scalars or Wigner spin 1 3-vectors. 

ii) the fourth one identifies the invariant mass of the isolated system as the Hamiltonian 
of the evolution in r when, with a gauge fixing, r is made to coincide with the Lorentz scalar 
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rest-frame time T, = Ps ■ ^s/ yPs = Vs ■ Xs/ \Jp1 of the decoupled external center of mass. 

In this description the standard manifestly covariant fields like the Klein-Gordon field 
0(2;^') are replaced by the new fields (f){T,a) = (j){z^{T, a)), which know the embedding and 
have the non-local information about the equal time hypersurfaces S(r) built-in. In the 
case of gauge theories, one can make a canonical reduction to a canonical basis of Dirac's 
observables in the radiation gauge (or Coulomb or generalized Coulomb; the literature is 
ambiguous about the terminology to be used), with the only universal breaking of manifest 
covariance connected with the external center of mass, since the relative motions are Wigner 
covariant. 

See Ref. for a complete review of the research program aiming to give a unified 
description of the four interactions in terms of Dirac-Bergmann's observables in the frame- 
work of the rest-frame instant form of dynamics, which is the classical background of the 
Tomonaga-Schwinger formulation of quantum field theory. 

The description of scalar (or spinning [P]) particles on arbitrary spacelike hypersurfaces 
requires the choice of the sign of the energy of the particle. This happens because the 
position of the particle on S(r) is identified by 3 numbers a = fj{T) and not by 4: this implies 
that the mass-shell first class constraints of the standard manifestly covariant approach have 
been solved and that one of the two disjoint branches of the mass spectrum has been chosen. 
In this way, one gets a different Lagrangian for each branch of the mass spectrum of an 
isolated system of particles (in the standard manifestly covariant description the branches 
are topologically disjoint for free particles): there is no possibility of crossing of the branches 
(the classical background of pair production) when interactions are present inside the isolated 
system (for instance charged particles plus the electromagnetic field) as happens in the 
manifestly covariant approach. While there is no problem in the coupling to magnetic 
fields of these particles with a definite sign of the energy, the minimal coupling in the 
Lagrangian will miss those couplings to the electric fields which are at the basis of the non- 
diagonalizability of both the Feshbach-Villars description of the Klein-Gordon field and of 
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the Dirac equation through the Foldy-Wouthuysen transformation (in the case of spinning 
particles these couphngs will have to be extracted from the iterative diagonalization of these 
theories and added non- minimally). However, this description of particles with, say, positive 
energy (whose quantization requires pseudodifferential operators |T^) seems suited for the 
description of the asymptotic Tomonaga-Schwinger states and will be used to introduce 
a notion of particle in a future quantization of classical fields on Wigner hyperplanes in 
the rest-frame instant form. These asymptotic states will replace the Fock ones of the 
standard manifestly covariant theory, which are the main source of problems in the theory 
of relativistic quantum bound states (the spurious solutions of the Bethe-Salpeter equation; 
see Ref. This framework should allow the introduction of bound states among the 



asymptotic states. 

Coming back to the isolated system formed by scalar charged particles of positive 
energy plus the electromagnetic field, we recall that in Ref. after canonical reduction 
to the radiation gauge, the final invariant mass of the reduced system is a function only of 
Dirac's observables (gauge invariant particles dressed with a Coulomb cloud and transverse 
radiation field) and contains: 

i) the kinetic energy of the radiation field; 

ii) the kinetic energy for the particles with minimal coupling to the radiation field; 

iii) the instantaneous action-at-a-distance Coulomb potential among the charges with the 
Coulomb self-energies regularized (the i ^ j rule, QiQj ^ 0) due to the Grassmann character 
of the electric charges Qi , i.e. = [at this semiclassical level we have Qi = eO*Oi] this 
does not imply the vanishing of the fine structure constant a = e^/Air ^ 1/137, since it 
gets contributions from QiQj; therefore, we are retaining effects of order a but not of higher 
order, because, as it will be shown, we do not have many-body forces]. 

Then, in Ref. |]12[ there was the evaluation of the retarded Lienard-Wiechert poten- 
tials of the charged particles in the radiation gauge in the rest-frame instant form. Since 
the Lienard-Wiechert potentials and fields are linearly dependent on the charges Qi, the 
semiclassical regularization Qf = eliminates the radiation coming from a single particle 
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(the electromagnetic energy-momentum tensor contains only terms in QiQj, i 7^ j) and the 
causal problems of the Abraham-Lorentz-Dirac equation of each charged particle [since the 
radiation reaction term has the coefficient Tq = 2Q^/3mc^ (tq is proportional to the time 
needed for light to travel across a classical electron radius), which vanishes if = 0, there 
are neither the Schott term nor the Larmor one but only the Lorentz forces produced by the 
other particles]. In Refs [|l3l one may find an extended discussion about this equation and in 



Ref. a recent review on its derivation and its causal problems (preacceleration, runaway 



solutions). See Ref. |jT5[| for modern attempts to extract the subset of causal solutions of this 



equation with the requirement of selecting only its solutions which admit a smooth limit 
for To ^ (the runaway solutions are singular in this limit) and to find an effective second 
order equation for this subset of solutions. 

Even if at the semiclassical level the "single charged particle" has no acausal behaviour, 
because, notwithstanding it produces a Grassmann- valued vector potential, it does not irra- 
diate, we can recover the asymptotic Larmor formula for a system of charges (considered as 
external sources of the electromagnetic field) due to the interference radiation from QiQj, 
i ^ j, terms (this result is in accord with macroscopic experimental facts). See Ref. ^ for 
the extension to spinning particles. 

Being in the radiation gauge, at each r the retarded Lienard-Wiechert potentials evalu- 



ated in Ref. contain also a non-local (in (j)term (coming from the transverse projector; 



see Eqs.( p.l3D , (|5.22|) in Section V): since this term involves all the points of S(r), the 



Lienard-Wiechert potential receives contributions from "all" the retarded times before r, 
namely from the whole past history of the particles. In the absence of incoming radia- 
tion one could put these retarded Lienard-Wiechert potentials inside the Lagrangian given 
in Eq.(74) of Ref. for the description of the isolated system in the rest-frame radiation 
gauge on the Wigner hyperplanes: one would get a Fokker-like action in the radiation gauge, 
replacing the standard Fokker- Tetrode one in the Lorentz gauge of the manifestly covariant 
description, and would have to face the problem of how to find a Hamiltonian description 
when there are integro-differential equations of motion with delay. The existing attempts 



are based on the idea of replacing retarded particle coordinates and velocities with instan- 



taneous coordinates and accelerations of all the orders (see for instance Refs. |T6|-pO|). In 
this way one replaces integro-defferential equations of motion with an infinite set of cou- 
pled differential ones. Since it is unknown how to formulate the Cauchy problem for these 



integro-differential equations (see Ref . ; exceptions are the 1-dimensional case or the 
time-asymmetric case [^,^), there have been complicated attempts to find conditions for 
extracting a set of effective second order differential equations from the infinite set (see for 



instance Refs. |l21| , |20| ). In Ref. there was an attempt to study the Dirac constraints 



originating from actions depending on accelerations of all orders, following previous attempt 



of Kerner [Id] of defining a Hamiltonian approach. See also the recent approach of Ref. [27 



Moreover, one should face a problem similar to the one raised in Ref. p8[, that only with 



symmetric Green's functions like ^{retarded + advanced) can the Fokker- Tetrode action 
corresponding to the Lorentz gauge give rise to a variational principle whose extremals are 
equivalent to the subspace of extremals of the original action defined by the symmetric 
Lienard-Wiechert solutions without incoming radiation, i.e. the adjunct Lienard-Wiechert 
fields (otherwise there are problems with the boundary terms) . This is compatible with the 
Feynman- Wheeler starting point for their theory of the absorbers (see Ref. |^ for the 



definition of radiation in this theory). A noncovariant justification of the results of Ref. ||28 
is given in Ref. [0; by ignoring the self-interactions and assuming a Lagrangian for two 
charged particles at equal times in which each particle interacts with the retarded Lienard- 
Wiechert potential of the other one, one obtains in the equations of motion Lorentz forces 
which correspond to ^{retarded + advanced) interactions, because in this Lagrangian the 
transition from retarded to ^{retarded + advanced) interactions is a total time derivative. 
However, self-reaction is ignored in these calculations and it is not clear how to arrive at a 
covariant formulation of these results. Let us remark that from the point of view of quantum 
field theory its regularization and renormalization require the use of the complex Feynman 
Green function (which does not vanish outside the lightcone; the solutions with the retarded 
Green function cannot be regularized at the distributional level): while its imaginary part is 
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connected to absorption in other channels, its real part is just the ^{retarded + advanced) 
Green function like in Feynman- Wheeler [ see Ref. ||31| for the extraction of a Fokker- Tetrode 
action with ^{retarded + advanced) kernel from the particle limit of QED (it does not work 
in QCD)]. 



However, both in the Dirac derivation ||32| of the Abraham-Lorentz-Dirac equation 



through the evaluation of the near zone self-field (with the same results obtainable with 
balance equations using the far zone fields; see Refs. and in the Feynman- Wheeler 

approach with the assumption of complete absorbers [|l3l the radiation is determined by 
the radiative Green function ^{retarded — advanced) [at the conceptual level this intro- 
duces the acausal advanced Green function and interpretational problems]. Indeed, the 
regularization of the self-energy divergence due to radiation reaction is done by rewriting 
retarded = ^{retarded + advanced) + ^{retarded — advanced) , by noting that the non ra- 
diative Coulomb piece of the fields (which does not influence the motion of the particle only 
giving a divergent electromagnetic contribution to the mass) is in ^{retarded + advanced) 
and by discarding this term as a regularization. 

However, till now all the calculations have been done in the Lorentz gauge and it is not 
clear whether the previous statements are "gauge invariant" . 

Since we now have the results of Refs. ||l],|l2| in the radiation gauge and a new type 



of regularization with the semiclassical approximation, it is interesting to revisit all these 
problems. 

The aim of this paper is to show that, starting from the rest-frame instant form descrip- 
tion of charged scalar particles plus the electromagnetic field with Grassmann-valued 
electric charges Qi, the semiclassical regularization = allows one to transform the sub- 
space of Lienard-Wiechert solutions (without or with incoming radiation) into a symplectic 
submanifold of the space of all solutions. This takes place since all the higher accelerations 
coming from an equal time development of the delay decouple, being of order on the 
solution of the particle equations of motion. As a consequence, at this semiclassical level the 
retarded, advanced and symmetric Lienard-Wiechert solutions coincide by using the equa- 



tions of motion, so that there is only one sector of semiclassical Lienard-Wiechert solutions 
(modulo the incoming radiation). The semiclassical Lienard-Wiechert potential and fields 
can be expressed as phase space functions and it is possible to eliminate the electromag- 
netic degrees of freedom by means of second class constraints added to the reduced phase 
space of the radiation gauge in the rest frame. Having gone to Dirac brackets with respect 
to these constraints, we get a reduced phase space containing only particles with mutual 
instantaneous action-at-a-distance interactions. We can find new canonical variables for the 
particles corresponding to a Darboux basis for these brackets. We can evaluate the final 
Hamiltonian, showing that besides the Coulomb potential there are vector potentials under 
the particle kinetic energy square roots (coming from the minimal coupling to the radiation 
field) and a scalar potential outside them (coming from the energy of the radiation field): 
due to Qf = one can extract the vector potentials from under the square roots and write 
a unique effective scalar potential added to the Coulomb one. This can be done both in the 
original (no longer) canonical variables and in the final canonical basis. The effective scalar 
potential is the complete Darwin potential: at the lowest order in we obtain the known 
form of the Darwin potential. 

We find that in the framework of Maxwell (not Feynman- Wheeler) theory, the semiclas- 
sical regularization = 0, QiQj 7^ extracts automatically the instantaneous action-at-a- 
distance potential hidden in the delay, which as mentioned above turns out to be the same in 
the retarded and advanced solutions, because the difference is in the terms which depend 
on the higher accelerations (in QFT these effects are hidden in the radiative corrections 
coming out from the regularization of the ultraviolet divergences, and this is possible only 
with the Feynman Green function). 

This means that at this semiclassical level, the elimination of the electromagnetic degrees 
of freedom produces a system of particles with instantaneous action at a distance given by 
the Coulomb and Darwin potentials. Since ^{retarded — advanced) — 0, all the effects now 
come from the regularized retarded = advanced = ^{retarded + advanced) solution and 
there is no mass renormalization. Even if the transverse projector implies contributions 
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from the whole past (or future) history of the particles, in the semiclassical approximation 
only the instantaneous act ion- at- a- distance effects on Svi/(t) survive. Each particle feels only 
the action of the other — 1 [thus giving us an effective Abraham-Lorentz-Dirac equation 
with no self-reaction and with the Lorentz forces of the other particles replaced by action- 



of radiation only as the effect of the other — 1 particles on the one chosen as a detector 



accord with the Larmor formula coming from the QiQj interference terms; (the Larmor 
formula gives zero at the semiclassical level due to the particle equations of motion). 

An important and new feature of this formalism is the possibility to find the final canon- 
ical variables for the particles after the introduction of the Dirac brackets: their use intro- 
duce new higher order contributions to the Darwin potential coming from the kinetic energy 
square roots. These contributions lead to a substantial cancellation with corresponding 
terms coming from what began as the electromagnetic energy integral. Our final generalized 
Darwin interaction naturally divides into two portions. The first portion has the same form 
either as the original lowest order correction derived originally by Darwin in the retarded 
case or as the lowest order 1/c^ well known form of the Darwin potential in the case of sym- 
metric []^{retarded + advanced) ] Lienard-Wiechert potentials but with the masses replaced 



a higher order correction, but it shows that we are using the correct relativistic kinematics 
without 1/c^ expansions). The second portion, a double infinite series, is, like that of the 
generalization of the 1/c^ Darwin potential, is also new. It is of higher order in 1/c^ than the 
more familiar first portion. Furthermore, our generalized Darwin interaction for A^ bodies 
is equal to the pairwise sum of two body pieces 

For the restricted case of two bodies considerable simplifications result when one evalu- 
ates the series in the center-of-mass rest frame. It then can be written in closed form. 

The Darwin potential we obtain can be regarded as the classical analogue of the full 
effects (complete transverse as well as longitudinal) of the single photon exchange in the 



at-a-distance interactions]. Like in the Feynman- Wheeler approach [^, we can now speak 



of radiation when it is far away from the other particles: 



equations of motion), one has 




and this is a new result 
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Bethe-Salpeter equation since it is the same order in the couphng constants. The effect 
of the semiclassical regularization = is to truncate out the classical analogue of the 
numerous higher order ladder and cross ladder diagrams. To the extent that the Darwin 
potential we obtain has a low order (1/c^) portion that agrees with the standard result, it 
would be expected to contribute correctly to the spectral results in a quantized formalism. 

For two particles there is the problem of the comparison of the semiclassical Lienard- 
Wiechert sector with the 2-particle models defined by two first class constraints with covari- 
ant instantaneous act ion- at- a- distance interactions (phenomenological approximations of the 
Bethe-Salpeter equation). Several authors beginning with [^^ -|36|| and continued by [^^HH 



have developed pairs of commuting generalized mass shell conditions, first class constraints 
with instantaneous potentials in the center-of-mass sytem, whose quantization gives coupled 
Klein-Gordon equations for two spin zero particles [see Refs. [p8|-p]l for similar equations 



for Dirac particles deriving from pairs of first class constraints for spinning particles]. Some 
of these models were generated as approximations to the Bethe-Salpeter equation, by reduc- 
ing it in a covariant instantaneous approximations to a 3-dimensional equation (with the 
elimination of the spurious abnormal sectors of relative energy excitations) of the Lippmann- 
Schwinger type and then to the equation of the quasipotential approach [see the bibliography 
of the quoted references], which Todorov |^5| reformulated as a pair of first class constraints 



at the classical level. In Ref. ||39| it is directly shown how the normal sectors of the Bethe- 
Salpeter equation are connected with the quantization of pairs of first class constraints with 
instantaneous (in general nonlocal, but approximable with local) potentials like in Todorov's 



examples. Ref. [^ shows how to derive the Todorov potential for the electromagnetic and 



world scalar case from Tetrode-Fokker-Feynman- Wheeler dynamics with scalar and vector 
potentials [this theory is connected with ^{retarded + advanced) solutions with no incoming 
radiation (adjunct Lienard-Wiechert fields) of Maxwell equations with particle currents in 
the Lorentz gauge]; besides the Coulomb potential, at the order one gets the standard 
Darwin potential (becoming the Breit one at the quantum level when spin is added), which 
is known to be phenomenologically correct. 
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What are the connections between the center-of-mass rest frame form of the two-body 
interaction Hamiltonian that we develope in this paper to all orders in and that obtained 
in the above references? The Darwin potential becomes a common overlap of the two 
approaches and thus an important testing ground for the approach we develope in this 
paper. Furthermore, when in future papers pseudoclassical spin is introduced (extending 
as in 1^ to interacting systems of particles and fields) the types of tests we perform in this 
paper will be relevant (note, however, there are difficulties with other categories of Darwin 
type of interactions brought on by the introduction of transverse spin-dependent electric 
field effects which unlike magnetic fields cannot be diagnonalized by a Foldy-Wouthuysen 
transformation [^|9|). 

In Section II we give a review of parametrized Minkowski theories on arbitrary spacelike 
hypersurfaces. 

In Section III we apply this formalism to the isolated system of charged scalar particles, 
with Grassmann-valued electric charges, plus the electromagnetic field and we arrive at its 
rest-frame instant form on the Wigner hyperplanes. We also make the canonical reduction 
to the radiation gauge. 

In Section IV we study the "internal" Poincare algebra and the "internal" center of mass 
on the Wigner hyperplanes and we derive the Hamiltonian and Lagrange equations of motion 
for fields and particles. Also the energy-momentum tensor is evaluated. 

In Section V we evaluate the Lienard-Wiechert potentials, we show that at the semi- 
classical level they depend only on particle coordinates and velocities (since at this level 
retarded = advanced) and we find their phase space expression. Then we eliminate the 
electromagnetic degrees of freedom by means of second class constraints which force them 
to coincide with the Lienard-Wiechert solution. We introduce the associated Dirac brackets 
and we find their canonical Darboux basis. 

In Section VI we derive the physical Hamiltonian and the effective Darwin potential to 
all orders in in terms of the old (noncanonical) and of the new (canonical) variables. In 
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the two-particle case we get a closed form of this potential using the rest-frame condition. 
Also the final form of the energy-momentum tensor is given. 

In the Conclusions there are some general considerations and hints for future develop- 
ments. 

Appendix A gives an explicit summation for the vector potential presented in Section V 
from the Lienard-Wiechert series in the rest-frame radiation gauge. 

In Appendix B we evaluate the field energy and momentum integrals used in Section VI 
when the electric and magnetic fields are expressed in terms of the Lienard-Wiechert series 
for the vector potential. 

We derive in Appendix C a general formula for a certain quantity, which is important 
for obtaining the closed form expression of the Darwin potential of Section VI for N=2 in 
the rest frame. 

In Appendix D we use a technique similar to that developed by Kerner ( and applied in 
p8| to obtain the Todorov quasipotential from the Wheeler-Feynman action) in the transfor- 
mation of the Lagrangian expression for the invariant mass to M. In this proof, (done to all 
orders) we must use Dirac brackets since we have used the Lienard-Wiechert constraints as 
a strong condition on the dynamical variables. This necessitates the explicit expression for 
the field momentum integrals developed in the earlier appendix from the Lienard-Wiechert 
solutions. 

In Appendix E we obtain a special solution of Hamilton's equations for the two-body 



problem that is analogous to Schild's solution for circular orbits [|2 
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II. PARAMETRIZED MINKOWSKI THEORY FOR N FREE SCALAR 

PARTICLES. 



In this Section we shall review the description of scalar free particles on arbitrary 
spacelike hypersurfaces j^, leaves of the foliation of Minkowski spacetime associated with 



one of its 3+1 splittings following the suggestions of Refs. ||4^,Q. The scalar parameter 
r labelling the hypersurfaces S(r) allows one to introduce a covariant concept of "equal 
time" , which will be useful in the description of an isolated system of interacting particles 
and fields. It will be shown that in these parametrized Minkowski theories there are first 
class constraints implying the independence of the description of the isolated system from 
the chosen 3+1 splitting. 

As said in the Introduction this requires the addition to the theory of an infinite number of 
new configuration variables z^(r, a) describing the spacelike hypersurfaces as an embedding 
of into Minkowski spacetime. We use the notation = (cr^ = t, a'"), i.e. A = (r, f) 
[the notation A = (r, r) will be used for the Wigner indices on the Wigner hyperplane, see 
Section III]. 

In the manifestly Lorentz covariant approach the worldlines of scalar particles are de- 
scribed by 4-vector coordinates a;f (xj), where the Xj's are affine parameters (often they are 
restricted to be the proper times of the particles). Even if one uses a unique affine param- 
eter T for all the particles, (r), there is the problem that the particles times x^{t) are 
gauge variables due to the presence of the first class mass-shell constraints pf — mf ^ 0. 
For each free particle the constraint manifold is the union of two disjoint submanifolds 



= ±ymf +Pi. The gauge nature of the x^'s is connected with: i) the arbitrariness in 
the choice of the center-of-mass time; ii) the arbitrariness in the choice of how to trigger 
the N particles (at equal times or with any conceivable mutual delay; this is the gauge 
freedom of relative times). Given the foliation of Minkowski spacetime with leaves S(r) we 
can give a covariant description of the particles at "equal times" (covariant zero relative 
times condition) by parametrizing the worldlines as a;f(r) = z^{rif{r)) = z'^{T,fji{T)), with 
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with the Lorentz scalar coordinates r]f{T) — (r, 77[(t)). Only 3 Lorentz scalar coordinates 
f/j(r) identify the intersection of the particle worldhne with E(t). This imphes that in this 
description there are no mass-shell constraints, namely that we are describing particles with 
a well defined sign of the energy: r)i — signp^ — ±1. In this paper we shall consider only 
positive energy particles, so that rji = 1 for every i. There will be a conjugate momentum 
Ki{T) for each particle and the standard momentum (r) will be a derived quantity which 
satisfies pf — mf. 

The metric induced on S(r) from the Minkowski metric r]'^" = (H ) is 

9AB{z{r, a)) := ri^,^z^^z% (2.1) 
where we have used the notation 



-l^-^-) = ^A--=^AZ^- (2-2) 



If we define the quantity z^{t, a) by means of 

44 = 4, (2.3) 

they satisfy 

9ab4''^ = Vxkz\zIz^z^ = T]^,. (2.4) 

Therefore, the ^^(r, a) are a set of vierbeins, with the z'^{t, a) the inverse vierbeins. 

Since we require > as a condition on the embedding z^ (r, a) , 2;^ is a time-like 
4- vector and the Zf{T, a)'s are spacehke 4- vectors tangent to S(t). 

The determinant of the metric is 

g ^ -det\\g^^\\ = {detz'^y. (2.5) 
The spatial part of the metric has an associated determinant which is defined by 

^ -det\\gfs\\; r = ^. (2.6) 
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We next define the inverse metric g by 



/^^BC = S^. (2.7) 

This imphes 

g'^^'gAB = 4 = /^447y^. = v'^v,., (2.8) 

so that 

= = + ^g"'"z'^z; + • (2.9) 

By definition of the element of an inverse matrix 

= (2.10) 

g g 

while the inverse of the spatial Qf^ is defined as 7''*, that is, 

Y'gst^sl. (2.11) 

To find g'^'^ use the fact that 

g^'grs + g^^grs^ST^o, (2.12) 

and therefore 

g^'grs = -g^^grs. (2.13) 

Multiplying by 7*" leaves us with 

= (2.14) 

Now using this consider 

g''gu + g''gtr = s!. (2.15) 

Multiply both sides by 7*", use the definition of 7*" and the above expression for g^'^ , and 
we obtain 
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= 7'" + -^r.^ri)7'V"- (2.16) 

9 

Thus in summary we have expressed the inverse metric in terms of the metric and the 
inverse of its spatial parts 

9'^ = -, 
9 

r2 

9^'' = ^7r.7^^ 

9 

grs^Y' + —9ru9rvl'''^Y'. (2-17) 

9 

Moreover, we have 

^^lu^AS _ fibu u A B _ Ab 



V'-'zX = = (2-18) 



The normal to E(t) at the point z^{t, a) is the Lorentz four vector 



nr, a) = f7^e^"''^^ia(T, a)zi^{T, a)zs,{T, a). (2.19) 



with the normalization 1^{t, a) — 1. By construction we have Z^(t, a)Zf{T, a) — 0. 

The evolution 4- vector z(^{t, a) can be decomposed on the mutually orthogonal four 

vectors and z'^: 



s 



z'^ir, a) = iV(T, a)P(T, a) + N'{t, a)4(r, a). (2.20) 

where N{T,a) is the lapse and N^{t,(j) the shift functions in the terminology of ADM 
general relativity. To determine the vector we use the orthogonality of 2;^^ and l^. That 
is multiplying both sides of the above equation by and using the expression for g^^ and 
multiplying by we obtain. 

N'^{r,a)^grr{r,a)Y''{r,a). (2.21) 

Using this expression for A'"'' we determine the scalar N by first multiplying the previous 
equation by z^r- Then we use the definition of and the determinant, multiply the result 
by g'^'^. Then use the definition of 7*"* , and we obtain N — T/y/j. Hence 



<(r,a) 



r 



V7 ^ 



(2.22) 



Substituting this expression for into 



(2.23) 



together with those for g^^ , we obtain after some algebra the following decomposition of 
the Minkowski metric 



= F{t, a)r{T, a) + Y'{t, a)z^{T, a)4(r, a). 



(2.24) 



Coming back to the scalar particles, the relation between their world line velocities in 
the two descriptions is 



(2.25) 



Noting that 



the standard action for N free scalar particles becomes 



(2.26) 



TV 



s=jdrY:[- 

with the Lagrangian density 



drL{T) = / dTd^a/:{T,a), 



(2.27) 



TV 



^r, a) = - m{r))mjgrr{r, a) + 2^, ,(r, a)r)t{T) + g ,,(r, a)7)[(r)?i|(r), (2.28) 



and the Lagrangian 



TV 



L(r) = -E^^^V^-(^'^^M) + 25-(r,r;,(T))77[(r) +^,,(r,77,(r))7)[(r)7)|(r). (2.29) 

i=l 

The above action is invariant under separate r and a reparametrization. This leads naturally 
to constraints. 



19 



The canonical momenta are determined from the dependence on the "velocity" z!^{t, a) 
associated with the hypersurface and the particle velocities tjUt). 



Kir [ T ) 



'9tt{t, a) + 2grf{T, ff^i^) + 9fs{r, a)f)i{T) ?)• (r) 

dL{T) 



(Jri-jT. //,(-)) + (J r.sjr. /7/(~))//f (~) on^ 

rrii — , _ _ , (2.30) 

ygrrir, f]i{T)) + 2gr f{r, f]i{T))f]l{T) + Q fs{T, ffi{T))rjf {T)rif (t) 



Using the above we have 



2 fs _ 2ri Y'igrf + gfiVi){9Ts + gsuVi)- 



m .2 



grr + 2grfr]l + gsfTilr]! 

{grr-Y^grfgrs)- (2.31) 



'grr + 2grrVl + gSfVM 

Use the following two forms 



=1 V^Tr(r, ffiir)) + 2grf{T, f]i{T))r)l + gsfir, fji{T))rilrjl 

Zr-l = -e^''^^ Zr^Zio.Z2pZ^^ = (2.32) 

together with the square root of the above relation and g'^^ = — '^grfY^^'^'^ 

g^'gn = g^^grA - g^^grr = i - -grr, (2.33) 

g 

to obtain 



Pl.l'' = E^'(^ - Vi{r))^m^ - Y'l^iri^is- (2.34) 



j=i 



Using finally 



u. 5^ (a - r]i{T))mi{z^Zrf, + z^z-s^^rjl) 



1 V^rr(T, fji{T)) + 2c/^^(t, r/i(T))7^[ + ^^^(t, fli{T))ril Tj? 

N 

Y,S%a-r)i{T))Kri, (2.35) 



i=l 
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we obtain the form of the four primary first class constraints Ti^ following from r and a 
reparametrization invariance: 

N . 

(r, a) = (r, a) - (r, a) ^ 5^ (a - r/i (r ) ) - (r, a) Ki^. (r ) (r ) - 

i=\ 
N 

- Zr^ir, a)Y'{r, a) S%a - r^(r))«,, ^ 0. (2.36) 

i=l 

Assuming the following Poisson brackets 

H{t),k,,{t)} = -6,,61. (2.37) 

one can show the first class nature of the above constraints IQ. Since these constraints are 
solved in terms of one of the independent momenta [p^(r, a)], their Poisson brackets are 
exactly zero: 

{n^iT,a),n,iT,a')} = 0. (2.38) 

These constraints imply that the description of the system is independent from the chosen 
3+1 splitting of Minkowski spacetime. 

The standard particle momenta are reconstructed as pf = {\Jmf — Y^i^ifi^is] f^i) and 
satisfy pf = rrif. 

In the next Section we shall add Grassmann-valued electric charges to the scalar particles, 
we shall find the new action and constraints, eventually arriving at the rest-frame instant 
form. 
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III. N CHARGED SCALAR PARTICLES AND THE ELECTROMAGNETIC 

FIELD. 



In this Section we will extend the formalism of the previous Section to the case of an 
isolated system of N charged scalar particles plus the electromagnetic field. By using Lorentz 
scalar electromagnetic potentials and field strengths, which employ the covariant "equal 
time" concept associated with the spacelike hypcrsurfaces S(t), we define the action for the 
combined field and particle system, that, like the free system, is separately invariant under 
T and a reparametrizations. As in the free particle case we obtain four primary first class 
constraints H^ir, a) pa at each point a on each of the space-like surfaces S(r), implying 
the independence from the chosen 3+1 splitting. Moreover, there are the two additional 
first class constraints describing the electromagnetic gauge invariance of the theory. By 
using a gauge fixing condition which restricts the S(t) 's to hypcrplanes Sj/(r) and by 
using Dirac brackets, the embedding variables z^{T,a) are reduced to only ten ones. The 
original constraints H^{r,a) are reduced to just 10 first class global constraints on 
each of the hyperplanes. We then specialize to hyperplanes orthogonal to the total timehke 
four-momentum of the system with 6 new gauge fixings depending on the standard Wigner 
boost for timelike Poincare orbits. These hyperplanes, defined by the system configuration, 
are called the Wigner hyperplanes Svk(t). After having defined the new Dirac brackets, we 
remain with a decoupled "external" canonical non-covariant center of mass, with Wigner- 
covariant particle and field degrees of freedom on the Wigner hyperplane and with only four 
first class contraints. One of these four constraints identifies the invariant mass of the system 
as the effective Hamiltonian, while the other three define the rest-frame condition Ti.p{T) ~ 
(vanishing of the total 3-momentum inside the Wigner hyperplane) for the combined 
system of particles and fields. Finally, we make the canonical reduction to eliminate the 
electromagnetic gauge degrees of freedom, placing the formalism (including the Poincare 
generators) in the Wigner covariant rest-frame radiation gauge. Also we give the energy- 
momentum tensor of the full isolated system. 
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A. The Action, The Constraints and the Canonical Reduction. 



Let us now review the isolated system of charged scalar particles plus the electromag- 
netic field following Ref. Just as on the hypersurface S(t) the positive energy particles 
are described by coordinates fjiir) such that xf(r) = z^ij^fjiij)), so the electric charge of 
each particle is described in a semiclassical way by means of a pair of complex conjugate 
Grassmann variables 6'j(r), 6*{t) ^ satisfying [J, = /* = 6*6i is the generator of the UemiX) 
group of particle i] 

= ef = 0, e^e* + e*ei = o, 

9i6j = 6j9i, 9i9* = 9*9i, 9*9* = 9*9*, i 7^ j- (3.1) 

As said in the Introduction, the formal quantization procedure sends 9*, 9i into the Clifford 
algebra describing a two-level Fermi oscillator 6|, bi, and each Grassmann- valued electric 
charge Qi = e9*9i goes into eblbi with eigenvalues ±e. 

The standard electromagnetic potential A^{z) and the field strength F^y{z) are replaced 
on S(t) by Lorentz-scalar variables A^(r, ct) and F^j^{T,a) respectively, defined by 

A^{t, a) = z\{t, a)A^[z{T, a)), 

F^5(r, a) = d^A^ir, a) - d^A^{T, a) = z^^{t, a)z^^{T, a)F^Mr, ^)). (3.2) 

The new potentials A^{T,a) have built-in the covariant concept of "equal time" through 
their implicit dependence on the embeddings z^{t,cx). 

With (i^S^ the surface element of S(r) we have the following volume element of 
Minkowski space-time 

c/S = z^drd^T.^ = drz^l^Vd^a = ^fgdrd^a. (3.3) 

The action now depends on the configuration variables z^{T,a),Aj^{T,a), ffi^r), 9i{T) 
and 9*{T),i = 1, .., A^, and consists of a "kinetic" piece for the complex Grassmann charges 
/ f [^i*('^)^«('^) ~ the same particle kinetic piece as in the previous Section, the 

23 



kinetic term / d'^z{—\F^^'^F^y) for the electromagnetic field, and the field-particle interaction 
term / QiA^(xi(r))xf (r) dr: 



S = J drd^aCir, a) = j drL^r), 
J d^aC{T, a), 



i ^ • • 

i=l 
N 



i=l 



+g,(r)(A.(r,a) + A,(r,a)77[(r))]- 



Q,{r)=eei{r)e,{r). (3.4) 

Since ^^-(r, a) transforms as a r-derivative the action is still invariant under separate 
r- and a-reparametrizations as in the free case. In addition it is invariant under the elec- 
tromagnetic local gauge transformations and under the odd global phase transformations 
59i y-^ iaOi , generated by the Jj's. The Qi = eli are the constants of motion associated 
with this last symmetry [-^Qiir) = 0, where '=' means evaluated on the solutions of the 
Euler- Lagrange equations; from now on we shall write Qi instead of Qiir)]. 

Since the semiclassical approximation Qf = regularizes the Coulomb self-energy, the 



criticism of Rohrlich |]13[ to this action principle [that the minimal coupling term and the 
electromagnetic field term are ill defined because they diverge on the worldlines of the 
particles] does not apply. 

The canonical momenta are [Ef. = Ffr and Bf = ^e^^fFgi {e^^f = e^^^) are the "electric" 
and "magnetic" fields respectively; for g^^ — > r/^^ one gets vr'' = —Ef = E^] 
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-2zf (r, a)^^-(T, a)/^(r, a)]F^^(T, a)Fc;5(T, B) 



dL 

OOrAT-{T, a) 

^'{r, a) = = _^M=7'""(r, - 9r,l'''F,,){T, a) = 

ddrAf[r, a) ^g(j^ ^) 



dL{T) _ 
d rilir) 

grfjr, fji{T)) + gfsjr, f)i{r))rjf{r) 

ygrrir, f]i{r)) + 2grf{T, f]i{T))f]l{T) + gfs{T, f]i{r))r]l{T)ril{T) 
+ QiAf{T, ffiir)), 



oei[T) 2 

d Oi (r) 2 
The following Poisson brackets are assumed 

{2;^(r, a), p,(r, a } = -r/f;(5=^(a - a ), 
{A4(r,a),7r^(r,a')} = r?|5^(a-a'), 

{^I('r),'«js(T)} = 
{^i(T),7rej(T)} = 

{^*(T),7r,.,(T)} = -5,,-. (3.6) 

The Grassmann momenta give rise to the second class constraints Trgj + f=:i 0, tt^i* , + 
^Oi ~ [{iiBi + ^9*,ng*j + |6'j} = —iSij]; ixei and Trg)*^ are then eliminated with the help of 
Dirac brackets 
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{A, BY = {A, B} - t[{A, 7le^ + '^e*}{7re^ , + '-9,, B} + {A, tt,. , + ^^,}{7r,, + '-9*, B}], 

(3.7) 

so that the remaining Grassmann variables have the fundamental Dirac brackets [which we 
will still denote {., .} for the sake of simplicity] 

{^,(r),^,(r)} = {^*(r),^;(r)} = 0, 

{9,{T),9*{T)} = -t6,j. (3.8) 
As in the free particle case of Section II, we obtain four primary constraints 
n,iT, a) = p,iT, a) - l,iT, a) [7;,(t, a) + J2 ^'i^ - Viir)) x 

i=l 



^Jml - Y^{t, a)[Kif{T) - QiAf{T, a)][Kis{T) - QiAs{T, a)] 



N 

(r, a)7^^(T, a){-Trs{T, a) + ^ 5\a - Ut))[k, , - Q^As{T, a)]} ^ 0, (3.9) 

i=l 



where 



2 a/7 2 
Trsir, a) = -F^f{T, a)7r*(r, a) = -e^ M7r*(r, a)B^(r, a) = 

^[n{T,a)xB{T,a)]„ (3.10) 

are the energy density and the Poynting vector respectively. We use the notation {nxB)s — 
{E X B)s because it is consistent with e^f^^n^By^ in the flat metric limit g^^ — > 77^^; in this 
limit -^1{E^ + B^). 

This form of the constraint displays both the tangential 2;^(t, a) and normal P(t, a) 
components of the momentum p'^(t, a) conjugate to the embedding variables z'^{t, a). 

Again, being solved in terms of the momenta p^(r, a),these constraints are first class 
with exactly zero Poisson brackets ({7i^(T, ct), 7i,^(T, ct')} = 0) and their existence implies 
once again that the description of the system is independent of the choice of foliation. 

Moreover, we have the (Lorentz scalar) primary constraints of the electromagnetic field 
connected with the gauge invariance of the action 
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n^{T,a)^0. (3.11) 

Prom these constraints we construct the Dirac Hamiltonian. But first we must construct 
the canonical Hamiltonian He- The canonical Hamiltonian is 

Hc^-Y. ^^f{r)f}l{T) + / rfV[7r^(T, a)drA^{T, d) - p^(r, a)<(r, a) - £(r, a)] = 

i=l ■' 

= y dV[9^(7r^(T, a)A^(r, a)) - A^(t, a)r(T, ^)] = - / (^VA^(t, a)r(T, a), (3.12) 

after the elimination of a surface term. 

Note that because of the r and a reparametrization invariance, He nearly vanishes, 
except for the portion involving 

AT 

r(T, a) = drn'ir, a) - ^ Q^S^ia - ff.ir)). (3.13) 

i=l 

Thus the Dirac Hamiltonian is (A'*(t, a) and A^(t, a) are Dirac multipliers) 

//ij = y dV[A'^(T, a):^^^(T, a) + A,(t, a)7r^(T, a) - A,(t, a)r(T, a)]. (3.14) 

The requirement that the primary constraints be r independent {{7r'^{T,a), Hd} 0, 
{H^{t, (j), ifz)} 0) leads only to the Gauss's law secondary constraint 

r(r,a)aiO. (3.15) 

Since the embedding variables z'^ (r, a) are the only configuration variables with Lorentz 
indices, the ten conserved generators of the Poincare transformations are: 

j^. ^ j^. ^ J ^3^(^Mp- _ ^-p/^)(r, a), (3.16) 
(the subscript s stands for hypersurface variable) . Prom the first of these we obtain 

{z^iT,a),p':}^-r\ (3.17) 
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We can restrict ourselves to foliations whose leaves are spacelike hyperplanes ^h{t) with 
constant timelike normal {drb^ = 0) , hy imposing the following gauge fixings [x^ is an 

arbitrary origin] 

C^{t, a) = z^{t, a) - <(r) - 6^(r)a^ ^ 0. (3.18) 

In this expression bf{T), f — 1,2,3 are three orthonormal vectors, such that the constant 
and future pointing normal to the hyperplane ^h{t) [6^ = (6^, bf) are orthonormal tetrads] 
is 

P(r, a)^F^b^,^ e^,^^bUr)bl{T)bl{T). (3.19) 

Using the definitions of the vierbeins we obtain from the above gauge fixing the simplifica- 
tions 

z^(T,a)^b'^{r), 
z^{T,a)^x^{T)+b>^{r)a', 

9fsir,a)^-5f,, Y'ir,a) ^ -S'', ^{T,a)^l, (3.20) 
as well as a natural decomposition of the Lorentz generators into orbital and spin portions 

Si^^ = 6^f(T) j d^aa'p^ir, a) - b',{T) j (faa'p^ir, a). (3.21) 

Here S^'^ is the spin part of the Lorentz generators. 

These gauge fixings have the following Poisson brackets with the primary constraint 

{e{T,a),n.{T,a)} = -r^l^6%a-a'). (3.22) 

Therefore, we get a continuum set of second class constraints. They can be ehminated by 
forming the Dirac brackets 
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{A BY = {A, B}- j [{A C^(r, a)}{H,{T, ^).B}- {A, H^(r, a)}{C(r, a), 5}]. (3.23) 
For example, one finds that 

{x^sM* = -V"'- (3.24) 

In this way the infinity of continuum hypersurface degrees of freedom z^(r, ct), p^(t, a), are 
reduced to 20 : i) 8 are a;^(r), p^; ii) 12 are the 6 independent pairs of canonical variables 
hidden in 6^ and S'^^ = J^^ — (x^p^ — x^p^) [they have the following brackets consistent 
with the orthonormality of the tetrads h\ 0: = 0, {^r^^il = V'^'^^a " ^''^^"a^ 

{^r, Sf} = C^r^Sf, where C^^^ = 6;6fr]''^ + 6i;6^sV"'" - S^S^V""" - ^^^sV"'^ are the 
structure constants of the Lorentz algebra]. 

It can be shown that the following 10 first class constraints survive at the level of the 
Dirac brackets [they are 10 combinations of the primary constraints whose gauge freedom 
is not fixed by the gauge fixings ( p.l8| )] 



N 



1=1 

+ yd'a[n' + B']{T,cr)] + 

N 

+ bUr)[T.[^if{r) - QMT,ffi{T))] + J d^a[7r x S],(r,a)] ^ 0, (3.25) 



and 



+ ljd'aa'[n^ + B^]ir, a) 

- [&^(r)6r(r) -6::(r)6^(r)][E^[(r)[<(r) -g.A^(r,r/,(r))] + 

i=l 

+ J rfVa"[7f X BY{T,a)] ^ 0. (3.26) 



i=l 



N 
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These constraints say that coincides with the total 4-momentum of the isolated system 
and that 5*^*^ is determined by its spin tensor. 

The Dirac Hamiltonian becomes 

Hd = \^{T)m{r) + VW^'^'W + /c?V[A,(r,a)7r-(r,a) - A,(r, a)r(r, a)], where, due 
to the associated Hamilton equations, the new Dirac multipliers have the following interpre- 
tation [§: >(r) = - x^(r), A^''(r) = -^^(r) = i[6^(r)65:(r) - 6^(r)6^(r)]. 

Restricting our considerations to configurations with > 0, we make a further canonical 
reduction to the special foliation whose hyperplanes are orthogonal to p^. These hyperplanes 
are intrinsically determined by the system itself and are called the Wigner hyperplanes 
Svy(r). They can be identified fl]] by requiring the gauge fixings &^(t) ~ Lf^^=A{Ps,Ps) for 

~ o 

the constraints Ti^'^^r) ^ 0, where L^t,{ps,Ps) is the standard Wigner boost for timelike 
Poincare orbits. This implies /'^ = ^ Ps/\/l^- 

The rest frame form of a timelike fourvector p'^ is P = rjy^^l; 0) = rj^°rj^Jp^, P"^ = p^, 
where t] = signp^. Since we restricted ourselves to positive energy particles, rji = +1, we 

o 

shall put 7] = 1. The standard Wigner boost transforming P ^ into p^ is 

L^ip^p) = t>:,{u{p)) = 

M , ^p'^k {p^ + p^){Pu + Pu) _ 

~ Vu ' 2 o 

/' p- P +p2 

uf'ip) + u^'{p)){u^{p) + u^{p)) 



< + 2u^'{p)u,{p) - 



1 + 



i2 



z/ = e^{u{p)) = uf^ip) = p^/^p 

y = r e^{u{p)) = i-urip); K - ^y^^)- (3-27) 

o o 

The inverse of L^^{p^P) is L^j^{p^p)^ the standard boost to the system rest frame, defined 

by 

L^(p,p) = L/(p,p) = L^(p,p)|p.__p^. (3.28) 

We also use these boosts to define the following vierbeins [the ei^{u{p)) 's are also called 
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polarization vectors; the indices r, s will be used for A =1,2,3 and o for A—O] 



4i'^ip)) = '^Aip), 



which satisfy 



4{u{p)yBi<p))=VB, 

ri^'' = e^^{u{p))ri^^eUu{p)) = u^{pX{p) -Y^e>^{u{p))e';{u{p)), 

r=l 

Vab = e|^(?7,(p))?7^^e^(M(p)), 

d d 
Pag^e'X{u{p))^Pa-^e^{u{p)) = 0. 

With the Wigner rotation corresponding to the Lorentz transformation A being 



o, M 



i?^(A,p) = [L{P,p)A-'L{Ap,P)] . 



1 

R'j{A,p) 



pi 



1\0 



((A-^)% - 1)P;3(A-1) 



(3.29) 



(3.30) 



(3.31) 



we have that the polarization vectors transform under the Poincare transformations (a. A) 
in the following way: 



6(f(«(Ap)) = (i?-VA^e^(«(p))). 



(3.32) 



These boosts can be used to obtain the further canonical reduction referred to above. 
This takes place in two steps: i) firstly, one boosts to the rest frame the variables 6^, 5'^'', 
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with the standard Wigner boost L^u{Ps-,Ps) for timehke Poincare orbits; ii) then, one adds 
the gauge-fixings 6^ — L^a{Ps-,Ps) ~ and goes to Dirac brackets. It can be shown [|l| that 
after this special gauge fixing the Lorentz scalar 3-indices f become Wigner spin 1 3-indices 
r. Therefore, we get a rest-frame instant form of dynamics with Wigner covariance. 

The Lorentz generators become J^'^ = x^p^ — x^p^ + S'^^ with S'^^ given in Eq.(59) of 
Ref. We now get 'H'^^ij) = 0, i.e. Si!'^ is forced to coincide with the spin tensor of the 
isolated system. 

If we define the rest-frame spin tensor 

N 

Y. 

i=l 



i=l 

- [V^vf - V^vt] [ I d'aa' [vr x fi],(r, a) + 

N 



i=l 



N 



i=l 
N 

-Y.Vt{T)[Kl{T)-Q,A^{T,f^,{T))]) + 
i=l 

+ j (a''[7f X BY{t, a) - a^[7r X BX{t, 



N 



i=l 

^ ' 73„„rr-;2 , d2i 



2 



d'aa'ir + B']{T,a). (3.33) 



it can be shown that the form of S^'^ implies that the rest-frame "external" Poincare gener- 
ators are UTI 



sir crs„s 

JT ■■ = ^IP's - (3.34) 

P°s+Vs\Pi 
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Only in this special gauge do we get the separation of a decoupled "external" canon- 
ical non-covariant center of mass described by the 4 pairs ^^(t), p^, of canonical vari- 
ables {{x'^,Ps}** = —T]^'^) identifying the Wigner hyperplane Svk(t) [see Eq.(59) of Ref. 

for the expression of 5;^(t) in terms of and of the spin tensor] and of the "inter- 
nal" Wigner- covariant canonical variables ffilr), K,i{T), A^(r, a), 7r^(r, ct) living inside the 
Wigner hyperplane and with the Dirac brackets coinciding with the original Poisson brackets 
( {v!{r),^,,{T)r = -h^l {A^{T,a),7,^{T,a')} =45\a - a')). 

As shown in Ref. one can replace x^iP^ with a new canonical basis Tg = Ps ■ 
^s/\It^s — Ps ■ ^s/\f[^s (it is the Lorentz-invariant rest frame time), £s = ^Jj^,Zs = 
^J^{xs — x°ps/p°),ks = u{ps) with Zs having the same covariance of the Newton- Wigner 
position operator under the little group 0(3) of the timelike Poincare orbits. The 3-position 
canonical variable z^/eg is the classical background of this operator and describes the decou- 
pled "external" 3-center of mass, whose 4-position is x^. 

In this special gauge there is no restriction on p^: the four velocity u^{ps) = Ps/\/]^ = ^'^ 
describes the orientation of the Wigner hyperplane with respect to an arbitrary Lorentz 
frame. 

We obtain the following form for the constraints 7Y^(r) ^ 0: 



n^^{r) = J d^aW^T, a) = - 

-«^b)(i|rfV[7f2 + fi2](r,a) + 



N 



i=l 

-e>^iuip)){ld'a[nxB]\T,a) + 

N 

+ J2l^iiT) - Q,l(r, ffi{T))Y) ^ 0. (3.35) 

i=l 

Their projections along the normal and the tangents to the Wigner hyperplane are 

n{T) = u^{ps)n,{T) = 



N 



J2^m^ + [K,iT)-Q,AiT,ff,iTW + 



i=l 
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+ -j£a[n' + B']iT,a))^0, (3.36) 

= Y.[^^ir) - Q^A{T,ff,{T))] + / rfVfvf X 5] (t, (x) ^ 0. (3.37) 
i=l 

The first one gives the mass spectrum of the isolated field plus particle system, while the 
other three say that the total 3-momentum of the charged particles plus fields vanishes 
inside the Wigner hyperplane Svi/(t). This condition is the rest-frame condition identifying 
the Wigner hyperplane as the rest frame of the isolated system. 

The Dirac Hamiltonian is now Hd = \^{t)1-L^{t) = X{t)1-1{t) — A(r) • 7Yp(r), with 
A(r) ^ -Xsf,{T)u''{ps), Xr{r) ^ -x^^(r)e^(M(p,)). 

The two additional electromagnetic constraints are 

7r^(r, a) ^ 0, 

r(r,CT)^0. (3.38) 

In the rest-frame instant form of dynamics on the Wigner hyperplanes they are Lorentz 
scalar constraints [Ar{T,(f) and 7r^(r, a) are Lorentz scalars, while A{t,(t) and ^^{T^a) are 
spin-1 Wigner 3- vectors]. 

We now eliminate the electromagnetic gauge degrees of freedom by decomposing the 
above spin-one Wigner 3-vector canonical field variables into their transverse and longitudi- 
nal components 

a - - 

A(r, a) = A^{T, a) - -d ■ A{T,a), 

d - 

^(r, a) = n^{T,a) - —d ■ 7f(r, a) ^ 

a N 

^ i=i 

with A = -d"^. 

We re-express everything in terms of the Dirac observables: i) A^{T,a), 7f_|_(r, a), 
{Al{T,a),fcl{T,a)} = -PJ_'{a)5^{a - a) [PJ_'ia) = S''' + ^] for the electromagnetic 
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field; ii) r/i(r), Kj(r) = /tj(r) + <9 ■ A{T,a) for the particles [they now become dressed 
with a Coulomb cloud]; iii) 9*{r), Oi^r), such that Qi = e6*9i = e9*9i. 

This is known as the Wigner-covariant rest-frame radiation gauge. Note that K,i{T) — 
QiA(T,f]i(T)) = Kii^r) — QiA±{T,f]i(T)). Using the Gauss law constraint d ■ n(r,(j) ^ 
SiQj5^(a — fji) with ^6^{a — ffi) = — l/(47r|(T — and integrating by parts we separate 
out the Coulomb portion of the rest frame energy from the field energy integral. A simi- 
lar procedure on the field momentum integral simplifies the rest frame condition (and the 
expression for the internal angular momentum in the next section). Thus we find that the 
reduced form of the 4 constraints is 

= - {E \lmj + {k,{T)-Q,A^{T,ff,{TW + 

i=l 

= e, - M ^ 0, 

Hp{T) = J^ikiir) + J c/V[^x X ^](r, a) ^ 0, (3.40) 

where eg = -^/pf and M is the invariant mass of the isolated system. The rest-frame condition 
does not depend any more on the interaction as it must be in an instant form of dynamics. 
This procedure not only extracts the Coulomb potential from field theory but also regularizes 
the Coulomb interaction due to the semiclassical property = 0. 

If we add the gauge fixing — r ^ we get A(r) = —1 and the Dirac Hamiltonian 
for the evolution in the rest-frame time is Ho = M — X{t) ■ 'Hp{T) [see the more accurate 
discussion after Eq.( |4.3| )]. 

The embedding corresponding to the Wigner hyperplanes in the gauge T, = r is 
2;^(r, (t) = x^(r) + e'^{u{ps))(^^ with the origin of the 3-coordinates given by x^(r) = 
x^(0) + u^{ps)t + e^(M(ps)) ^ dr' Wt') since x^(r) = u^{ps) + e'^{u{ps))Xri'^) [instead for 
the "external" center of mass we have ^^(t) = M^(ps)]. The final canonical variables are: i) 
Zg, kg (the decoupled "external" 3-center of mass); ii) r/j(r), t^iij) (the particle variables); 
iii) A_L(r, (t), vf_L (r,(T) (the transverse radiation field). They are still restricted by the three 
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rest-frame conditions 'Hpij) ^ 0. 



B. The Energy-Momentum Tensor. 



The Euler-Lagrange equations from the action ( |3.4j) ) are 

drjl dr drjl ' 
where we introduced the total energy-momentum tensor [rjfir) = (l;?7j(r))] 

-hn^-m) 'l-^- '1 , + 

i=l 



+ [F-^^f/ + ^/^F^^F^6](r, a). (3.42) 

When dj^l^z^ (r, cr) = as happens on the Wigner hyperplanes in the gauge — r ^ 0, 
A(r) = 0, we get the conservation of the energy-momentum tensor T^^{t, a), i.e. d^T^^ = 
0. Otherwise there is compensation coming from the dynamics of the hypersurface. 

On the Wigner hyperplanes the energy-momentum tensor becomes 







N 
















i=l 








N 




rj-irr 






[a — 






i=l 








N 
















i=l 








rl 





[<(r) - Q,A^{T,Ur))M{r) - Q.A^irMr))] _ 

-5"^[7r" + B^] - [tt'tt^ + B'-B']] (r, a). (3.43) 

Finally, after the canonical reduction, which eliminates the electromagnetic gauge degrees 
of freedom and the choice of the gauge Tg = r, we get 
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N 



i=l 



1 ^ f) 



2 

AT 



i=l 

N ^ 



i=l 



N 



[k\{T) - Q,Al{T,fj,{T))\[kl{T) - Q,il(T,77;(T))] 



1 ^8 2 



1=1 

[(^± + E Qri::^\^ - + E Q^i::^^^ - Vi{r))y + 



i=l 



=1 



+ 



(3.44) 
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IV. INTERNAL POINCARE ALGEBRA AND EQUATIONS OF MOTION FOR 
THE ELECTROMAGNETIC FIELD AND THE N CHARGED PARTICLES. 

In this Section we first build a realization of the Poincare algebra inside the Wigner 
hyperplane using results of the previous Section. Then, by identifying the Lorentz scalar 
rest-frame time with the invariant time r labeling the hypersurfaces S(r), we arrive at 
the Dirac Hamiltonian Hd = M — A(r) ■ T-Cp{T), in which the only gauge freedom left is the 
one associated with the rest-frame condition. We then obtain the Hamilton and Lagrange 
equations for fields and particles. Then we describe how to find the canonical "internal" 
center of mass for fields and particles on the Wigner hyperplane. The natural gauge 
fixing to the rest-frame conditions 7ip(T) ^ are g+ ~ 0: they imply A(r) =0 and the 
decoupling of the "internal" center of mass from the "internal" relative motions. In this 
way only the "external" decoupled 3-center of mass Zg remains (the Newton- Wigner-like 
3-position which replaces the 4-center of mass in the gauge Tg = t). A property of the 
particle accelerations of any order, which will be needed in the next Section, is derived. 



A. Internal Poincare Algebra 

In the rest-frame instant form of the dynamics there is another realization of the Poincare 



algebra besides the "external" one given in Eq.( 3.34 ). This is the "internal" realization built 
in terms of the variables living inside each Wigner hyperplane. The associated generators of 
this internal Poincare group are given by (for positive energies in the Wigner-covariant rest 
frame radiation gauge) 



N 
i=l 



Viint) = np = Kir) + J X B]{t, a)^0, 
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1=1 

N 



i=l 

N 1..N -, o 

1=1 j^i ^fi] '^^j 

+ Q,Jd'anl{r,a)c{a- rf,(r))] - i | (#1 + ^')(r, a), (4.1) 

in which k,+ {t) = Y^iKiij) and c{f]i — f]j) := — l/(47r|r/j — f]i\). The latter two Lorentz 
generators are determined as the components of the spin tensor Sf^ defined in Eq.( p.33D 
inside each Wigner hyperplane. 
The Dirac Hamiltonian is 

Hd = X{r)n{T) - A(r) • 7?p(r). (4.2) 

As aheady said, if we add the gauge-fixing 



implying that the Lorentz scalar parameter r labelling the leaves of the foliation of Minkowski 
spacetime with Wigner hyperplanes coincides with the rest-frame time of the decoupled 
point particle clock (the "external" center of mass) x^, its conservation in r will imply A(r) = 
— 1 so that, after taking the Dirac brackets associated with the second class constraints 
eg — M ^ and — r ^ (this eliminates and eg), the final Dirac Hamiltonian in this 
gauge would be Hd = — A(r) ■ Hp^r). However, if we wish to reintroduce the evolution in 
r = T, in this frozen phase space [containing the canonical variables Zg, kg, fji^r), n-iir), 
^±(t, c?), 7fj_(r, a)] we must use the Hamiltonian 

HD = M-X{T)-np{T), (4.4) 

because M = VJ^^£~^ is the invariant mass and the "internal" energy generator of the isolated 
system [it is like with the frozen Hamilton- Jacobi theory, in which the time evolution can 
be reintroduced by using the energy generator of the Poincare group as Hamiltonian]. 
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The only remaining first class constraints are the rest-frame conditions TCp^r) ^ 0. The 
Dirac multipliers A(r) describe the remaining gauge freedom on the location of the "internal" 
center of mass on the Wigner hyperplanes. In the next Subsection we will study the natural 
gauge fixings for these first class constraints. After this final canonical reduction the isolated 
system will be described by the decoupled "external" 3-center-of-mass canonical variables 
Zg, kg and by relative Wigner-covariant degrees of freedom on the Wigner hyperplane, with 
the invariant mass M as the Hamiltonian for the evolution in r = T,. 



B. The equations of motion for particles and fields. 



The Hamilton-Dirac equations associated to the previous Hamiltonian are 

Vm2 + (^,(r)-Q,A^(r,77,(r)))2 
^(\^ QiQk{ffi{r) - rffc(r)) 

"'^"^ fe4vrhl(r)-r/.(r)|3 + 

k+{T) + / c/V[^x X B]{t, a) ^ 0. (4.5) 



in which = means evaluated on the equations of motion. 
The Hamilton-Dirac equations for the fields are 

A±r{r, a) = -7r_Lr(r, a) - [A(r) ■ d]Aj_r{r, a), 
^l(r, a) ^ Ail(r, a) - [A(r) ■ 9]7fl(r, a) + 

-J2Q,Pl^{amT)S'ia-v,iT)). (4.6) 



The associated Lagrangian, obtained by means of an inverse Legendre transformation [|12 
is 



LRir) = ffii^T) ■ Ki(r) - J dVv4j_(r, a) ■ 7r±{T, a) - Hr{t) 



N 

^ [ - m^l - (^.(r) + A(r))2 + Q,%{t) + A(r)] ■ ^(r, r/,(r))] + 

i=l 
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1 X ^ QiQj 

2 ^. 471 \ Ur) - m(T) \ 



Here A(r) is now interpreted as a non-linear Lagrange multiplier needed to get the rest-frame 

. d OLr ^ dLR. OLr 



conditions Hp = V(int) ~ 0. Its Euler-Lagrange equations^|^^ = = yield 



d , f/j(r) + A(r) ^ > / ^ , ... o 



■■r 



-i^(r,a)-^{[A(r).9]il(r,a)} = 



= Ai:L(r, ^) + [A(r) ■ 9]{i^(r, a) + [A(r) • d]Al{r, a)} - 

N 

- Y.QiP±{m^{r) + X\r)]5\a-UT)\ (4.8) 

i=l 

and [these are the Lagrangian rest-frame conditions] 
.=1 V^l-(^(r)+A(r))2 

+ / + [A(r) ■ 9]il)](r,a) ^ 0. (4.9) 
The Lagrangian expression for the conserved invariant mass M = VJ^^^^ is 



j-i \ ^ \ ^ QiQj 



^=1 - iUr) + A(r))2 4vr | //^(r) - 7S(r) 



\ - 2 ~ 2 

+ / c^V-[E^ + B ]{T,a) = const. (4.10) 



Eq.( [4.8|) may be rewritten as 



i_^^.^W + A(r^^ 
Vl + (^^.W+A(r))2 

47r I r/,(r) -r/fe(r) |3 
+ Qa^±(r,r/,(r)) + (r/,(r) + A(r)) x%,r/,(r))], (4.11) 
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where the notation = —Aj_ — [A(r) • d]A''[_ = tt\_ has been introduced. 

Eqs.( [4.1ID and ( ^.81 ) are the rest-frame analogues of the usual equations for charged 



particles in an external electromagnetic field and of the electromagnetic field with external 
particle sources in which both particles and electromagnetic field are dynamical. Eq. ( |4.9|) 
defines the rest frame by using the total (Wigner spin 1) 3-momentum of the isolated sys- 
tem formed by the particles plus the electromagnetic field. Eq. ([4.10| ) gives the constant 
invariant mass of the isolated system: the electromagnetic self-energy of the particles has 
been regularized by the Grassmann-valued electric charges [Qf = 0] so that the invariant 
mass is finite. 

C. The "internal" center of mass and the last gauge fixing. 

The rest-frame conditions Hp = V(int) ~ show that there are still 3 gauge degrees 
of freedom among the reduced canonical variables ffiir), K,i{T), A±{T,a), 7T±{T,a) on each 
Wigner hyperplane. They correspond to our freedom in the choice of the point of the Wigner 
hyperplane that locates the "internal" 3-center of mass g+ of the isolated system. After the 
gauge fixing g+ only Wigner-covariant relative variables are left on the Wigner hyper- 
plane and there is no double counting of the center of mass [only the decoupled canonical 
non-covariant "external" one Zg, kg is left]. 

In Refs. PIJT^ there was a naive choice ff^ = Y^iLi Vi of the Wigner spin 1 3- vector 



conjugate to Hp = V^mt)- Then, after realizing that ^ does not imply A(r) = 0, 
in Ref. |^ a different choice (f+ was made by utilizing the group-theoretical results of Ref. 
1^: now the time constancy of the gauge fixings ~ implies A(r) = 0. Moreover, the 
nonrelativistic limit of g'^ is now the unique nonrelativistic center of mass. 

In this gauge we get the simplest description of the dynamics on the Wigner hyperplanes: 
a = Q implies that the "internal" center of mass is put at the origin x^(r) of the 

coordinates [z^^ij.a) = x^ij) + e^('u(ps))cr^ with x^ij) = x^(0) + u'^{ps)t] of the Wigner 
hyperplane. In this gauge the origin acquires the property i^(r) = u^{ps) and becomes also 
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the "external" Fokker-Pryce center of inertia of the isolated system |^ [in a future paper 
there will be a more detailed analysis of these problems]. 

In order to find g+ one must take advantage of the "internal" realization of the Poincare 
algebra inside the Wigner hyperplane. Ref. implies the following definition of the canon- 
ical "internal" 3-center of mass 



^{int) ■ 'P{int)'P(int) 



nut) ^ 

Imposing the rest-frame condition, it is seen that qj^ weakly coincides with the non- 
canonical "internal" M0ller center of energy In that same limit it is also equal to the 
"internal" Fokker-Pryce center of inertia defined by 

where 

S[int) = J{int) -q+X V{vnt) ~ 5^^- (4.14) 

With the gauge fixing condition ^ and with Ts = r one finds the following expression 
for the origin x^(r) of the coordinates on the Wigner hyperplane [a;^(0) is arbitrary] 

xf+^^{Ts) = a;^(0) + n^(p.)T,. (4.15) 

It can be shown that this coincides with the covariant noncanonical "external" Fokker-Pryce 
center of inertia y^(r). However, it is different from both the "external" center of mass 
x^(r) and the "external" center of energy of M0ller i?^(r). 
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Since = {g+,M — A(r) ■ Hp} = — A(r) ^ 0, there is no gauge freedom left and 

we could eliminate the variables g+, 'P{int) = T^p and look for a canonical basis of (Dirac 
observable) relative variables on the Wigner hyperplane. 

Instead of doing that [see Ref. [§]], in this paper we will go on to work with all the 
variables r/i(r), ^^(r), Aj_{T,cr), 7f^(r, a), but we shall restrict their equations of motion to 
the gauge A(r) = without explicitly introducing g+ ^ 0. 

The equations of motion for the particles and for the electromagnetic field then become 

d-^'^'^V^) ^.4vr|,l(r)-rf,(r)|3 + 

+ Q^[Mdr.m{r))+Ur) x %,rf.(r))], (4.16) 



□il(r, a) = l^(r, a) + Ail(r, a) = Jl{r, a) 

N 



i=l 



1=1 ^ 

TV 



i=l 

+ [d'^' ''^T_Z,^'h %a'-ff,{r)l (4.17) 



a — a' 



with 



n'^{a - a') = 5" - 3(a" - (x"'){a' - a")/ {a - <?') • (4-18) 



We point out that defining /3j(r) = fj^ij-) = ^^^-^ = ^ [t = ct, r]^{t) = fji^r)] even 
if we use everywhere c = 1, we have momentarily reintroduced it] and fi^^^ = d^^i/dr^ and 
writing the particle equations of motion as (no sum over i) 

^im,^M^) = + ft-C^I - Q.F., (4.19) 

we obtain 
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so that 




(4.21) 



Thus in general we will have for every h>l 



(^t^ — QiGi, 



(4.22) 



so that using the Grassmann property of the charges 



0, 



h>l. 



(4.23) 



This will lead to important simplifications later allowing us to drop acceleration dependent 
terms in the force. 

Due to the projector PJ_^{a) required by the rest- frame radiation gauge, the sources of 
the transverse (Wigner spin 1) vector potential becomes non - local and one has a system of 
integrodifferential equations (like with the equations generated by Fokker- Tetrode actions) 
with the open problem of how to define an initial value problem. 

The Lagrangian equations identifying the rest frame become 



N 




r 



(4.24) 
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V. ELECTROMAGNETIC LIENARD-WIECHERT POTENTIALS. 



In this Section we will study the Lienard-Wiechert solutions of the previous radiation 
gauge field equations in the gauge A(r) = in absence of incoming radiation by using the 
results of Ref. ||T^. We shall study the retarded, advanced and ^{retarded + advanced) 



Lienard-Wiechert potentials. Using the Smart- Winter [56] expansion for retarded and ad- 
vanced time dependence, we obtain an infinite series form of the retarded and advanced 
Lienard-Wiechert potentials depending on instantaneous accelerations of every order. It will 
be shown that the results of the previous Section imply that on the solutions of the particle 
equations of motion the higher accelerations decouple due to the semiclassical regulariza- 
tion = 0. We show that this implies that the ^{retarded — advanced) Lienard-Wiechert 
potential vanishes at the semiclassical level and that there is only one semiclassical Lienard- 
Wiechert potential: retarded = advanced = ^{retarded + advanced). This allows us to 
re-express the semiclassical Lienard-Wiechert potential in terms of particle canonical coor- 
dinates and momenta [the same can be done for the Lienard-Wiechert electric field, as it will 
be shown in the next Section]. Therefore, we get a Hamiltonian description of the Lienard- 
Wiechert semiclassical solution: this sector of solutions can be identified as the symplectic 
submanifold of the space of solutions of the electromagnetic field equations determined by 
two pairs of second class constraints, which force the electromagnetic field to coincide with 
the semiclassical Lienard-Wiechert one.. After having gone to Dirac brackets with respect 
to them, we get a reduced phase space with only particles and we find a canonical basis ff^, 
Ki for these brackets for arbitrary A^. In the new variables the rest-frame condition becomes 
J2iLi ~ O5 expected in an instant form of dynamics. 
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A. Grassmann truncated form of the advanced and retarded Lienard Wiechert 

Solutions. 



Here we develop the Grassmann truncated forms for the Lienard- Wiechert vector 
potential [see the next Section for the transverse electric and magnetic fields]. The 
^{retarded + advanced) solutions are given (for A(t) = 0) by [for the sake of notational 
simplicity we will use the notation ^^(t), A±{T,a), 7T±{T,a), instead of Ki{T), A±{T,a), 
'K± (r, a) , from now on] 

1 



5[{t - nf -{a- a^m{n)5'{a^ - fj,{n)) = 
= Pr(a) E ^ y dt,5[{t - t,f --{a- Uct,)m{ct,) 

N 

■■ -Y.Qi^^'^Siir.Bl (5.1) 



j=i 



in which we have put r = ct, Piij) = fj^ir) = ^— ^ and A±^ = A^ret (^±- = Aj_adv) 
for the retarded (advanced) solution. The equation for ti is c^{t — tiY — {a — f]i{cti)Y with 
the two solutions being 

11 T 

tj+(r, a) = -Tj+(T, a)^t- -rj+(rj+(T, a),a)^-- Ti+{T, a), 
c c c 

11 r 
ti_{T, a) = -Ti-{T, a) = t + -ri-{Ti-{T, a),a) = - + Ti_{T, a), (5.2) 
c c c 

for the retarded and for the advanced case respectively. The hght cone delta function is 

s[{t - nf -{a- ff,{n)r] = i5[(t - t,r - i(a - ff,{ct,)r] = 

^ S[Ti_-Ti+{T, a)] ^ S[ti. - r,_(r, a)] 

2|r - n - A(ri) ■ {a - Un))] 2|r - n - A(ri) ■ {a - ff,{n))\' 

The relative space location between the field point and the retarded or advanced particle 

position is 

- fji{ri±{T, a)) = ri±{Ti±{T, a), a) = ri±{Ti±{T, a), a)fi±{Ti±{T, a), a), (5.4) 
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and its length is related to the time interval by 



ri±{Ti±{T, a),a) = \a- f]i{Ti±{T, = cTi±{T, a) = \t - ri±(r, a) 



T - ri±(r, a) = ±cTi±(r, a) = ±ri±(r,±(r, a), a) 



The effective spatial interval is defined by 



— * 

pi±{.n±{,'T, = ri±{Ti±{T,a),a)[lT l3i{,ri±{T,a)) ■ fj±(ri±(r, a), a)]. 



In terms of these variables, the retarded, advanced and time symmetric solutions a: 



7i 47r 



N 



pi±[ri±[r, a), a) 

N 



Q^ 



i=l 



i=l 



TT 



pi+{Ti+{T,a),a) pi-{ri-{T,a},a) 



We use the Smart- Wintner expansion P6| , pr8| , p!9 



fi-Ti±) = f{T =F cTi±{Ti±{{T, a), a))) = /(r - [±ri±(ri±(r, a), a)]) 



fir) + E - 

k=l 

k=0 



k df{ 



r 



dr 



k\ dr'' 



(±r,(r,a))'=[lTA(r)-f,(r,a)]/(T) 



where 



fj(r, a) = ri(r, a)fj(r, a) = a - f]i{T) = r;±(r±(r, a), a)\r,^^r,a)-- 



and where the last line in Eq. (^.8|) is identical to the previous one since '^'''^^''^^ 



(It 



r^r, a). 

Hence we get 

AT 



oo l^^Y d^ 
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In order to evaluate the above derivatives we need the Leibnitz formula for the kth derivative 
of the product f{T)g{T) 



Thus we get 



fc=0 



k\ dr^^ ' A;! ^Qm\{k - m)\ dr'^ dr''- 



^nt^ m\(k-m)\ dr"" dr''-'^ ~ ^nt^n rn\h\ dr^ dr"^ ' ^' ^ 

m=U k=m ^ ' m=U h=0 



Using the notation fi'^^^^ — we obtain the following expression for the vector potential 

^) = E |:n^(^) E ^/?P(r)0.±,.(r, a), (5.13) 

j=l /i=0 



in which 



In order to display the result of the evaluation of the derivative we use the formula 
(with the summations restricted so that = Z^r = m ) to obtain 



W-,-)=E^EE 



(5.16) 



art V ^' / \ ^' 

In this expression the symbol o represents a scalar product between the tensors to the left 
and to the right with the summationX^r Ur — n indicating how the indices would be matched. 
Changing the m summation index to A; = m — n we obtain 

n=0fe=0"in2.. ni\n2\.- 

(5.17) 



<^i±Ar. ^) = E E E — 



d-rr'^''-\r,d) ^ A(r) A^^^(r) 
d?y ° 1! 2! 
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(In this latter summation J^r'^r = n, X]r™r = n + k.) 

Now we can take advantage of the Grassmann charges to significantly simplify the above 
multi-summations. As we have seen above, with a semiclassical Qi there are no accelerations 
on shell {QSf^^ = 0) in the equations of motion of the particle H\ since both the Coulomb 
potential and the Lienard-Wiechert Lorentz force on particle 'i' produced by the other 
particles, i.e. Qi[E±{T,f]i{T)) + /3j(r) x 5(r, f/i(r))], are proportional to Qi. Therefore, 
the full set of Hamilton equations ( [4.5|) , ([4.6| ) for both fields and particles imply that at 
the semiclassical level we have a natural "order reduction" of the final particle equation of 
motion in the Lienard-Wiechert sector [only second order differential equations]. 

One effect of this truncation is the elimination of multi-particle forces; all the interactions 
will be pairwise, in both the Lagrangian and Hamiltonian formalisms. This was to be 
expected since the rest-frame instant form is an equal-time description of the particle 
system: (acceleration- independent) 3-body,.. N-body forces appear as soon as we go to a 
description with no concept of equal time, like in the standard approach with first class 
constraints ^ . 



Thus the only contributing indices are n2 = = .. = 0, ui = n and our expression for 
the transverse vector potentials simplify to 

^ O 

Al^ir,^) = E|jPr(^)A^(r)0.±,o(r,cr) = 

- E^n(-)AWE^^ ^ — [-^ 

Alsir^a) = 1(A1^ + Al_)(r,a). (5.18) 
Since = we see that for odd n=2m+l we get 

' ^)''" = ™T(^tr = 0, (5.19) 



and this implies the equality of the retarded, advanced and symmetric Lienard-Wiechert 
potentials on-shell 

A^^s{r,^) = Al^{r,a)= E gpr(a)A^(r) 
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E 



'o(2m)! 



A(r) ■ 



Q2ra \ 



,2to-1. 



T, a). 



(5.20) 



Therefore, at the semiclassical level there is only one Lienard-Wiechert sector with a uniquely 
determined standard action-at-a-distance interaction. 

We use a tensor notation to write the transverse symmetric vector potential above as 



(5.21) 



Using the definition of the Coulomb projection operator 



dahdak \a' - a 



(5.22) 



— * 

and compactifying the notation still further we obtain [Vo- = d/da] 



N 



Q 



i=l m=0 l^^^'J- 

47r J \a — (T 



|2m-l 



(5.23) 



Integration by parts and changing from ^ to ^ and translation gives 



N 



Q 



1 



a±s{t, ^) = e ? e %{T){ffAr) ■ y.fn 1^ - m 



\2m-\ 



(5.24) 



The integral above is finite, and thus we can view it as the A — > oo limit of an integral with 
a cutoff A and take the derivatives out. The integral is thus of the form 

1 



/2m-l 



mm 



(5.25) 



and 



1 r _/2m— 1 

An J A \a' — {a — ffi 



2io J-i 



dz 



a' + {a — ffiy — 2a'\(T — ffilz 

r da'a""^^' , ^ , J a'' + {a - fj.f -2a\a- ff, 
Jo a'\a — rjil * 

1^"^ ^1 / da'a''^"'{\a' - \a - f]i\ \ - \a' + \a - f]i\\) 
\(J — rjil Jo 



1 /-A 

2 



^2m+l 



- Vi 



2m+l 



2m + 1 (2m + l)(2m + 2)' 



(5.26) 



51 



Note that the A cutoff will get killed by the a derivatives. Thus, we obtain 

N oo -| 



(2m + 2)! 

TV 



V.(^.(r)-V.)2™+i|a-r/, 



r 



■.=J2Q^^±s^i^-V^ir),f]Ar)). (5.27) 

i=l 

Using the first half of particle Hamilton equations ( |4.5| ) [with A(r) = 0] in the form ffi = 



K,i/ymf + K,f + 0{Qi), we can, as shown in Appendix A, arrive at the following closed form 
of the vector potential [ffi = ffi{T), Ri = Ki(r)] 



N 



i=l 



1 



-Ki ■ (I 



ATx\a — f]i 
{a - f]i){a - ffi) 



mf + (Hi ■ 



|o" - m 



)( 



mj + Ki^ 



mi + [Ki 



l)x 



mi + Ki 



K; I Ki 



(5.28) 



B. Lienard-Wiechert Second-Class Constraints, their Dirac Brackets and the New 

Canonical Variables. 



Thus far we have the reduced phase space of charged particles plus the transverse 
electromagnetic field. This is a well defined isolated system with a global Darboux basis [ 
r/j, Kj, y4_L(r, a), 7f^(r, a)] and a well defined physical Hamiltonian, the invariant mass M = 
'^(int)- possible configurations of motion take place in this reduced phase space. The 
space of solutions of Hamilton's equations is a symplectic space in that there is a definition 
of Poisson brackets on the space of solutions. The question arises whether one can select 
a subset of solutions of the equations of motion which is still a symplectic manifold: an 
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arbitrarily chosen set of solutions will not form a symplectic manifold. The method we 
propose here is to add by hand a set of second class constraints "compatible with the 
equations of motion" which amounts to the selection of a symplectic submanifold of the 
symplectic manifold of solutions. 

— * 

The above Grassmann truncated semiclassical Lienard Wiechert solution A±s for the 
vector potential with 7r_L5 = E±s = —^^±5 for the canonical conjugate field momentum 
[see Eq. (|6.2| ) in the next Section] and provide us such a set of second class constraints by 
way of 

Xi{t, a) = A±{t, a)-Y^ Qi^±Si{^ - Viir), Ki(r)) ^ 0, 

i=l 
N 

X2{t,B) = n_i{T,a) -^QiTiLSii^ - r]ii^),f^iir)) ^ 0. (5.29) 

1=1 

These constraints allow us to eliminate the canonical degrees of freedom of the radiation 
field and to get the symmetric Lienard-Wiechert reduced phase space, in which there are 
only particle degrees of freedom. This has an immediate and important consequence: the 
independent variables f/j, will no longer be canonical when one imposes these constraints 
by way of modified Dirac brackets. 

Now in order to compute the effects of these constraints we must use them in the con- 
struction of Dirac brackets. This requires that we compute the 6x6 matrix of brackets 



{Xi,Xi} {Xi,X2} 

{X2,Xl} {X2,X2}J' 



(5.30) 



It turns out that this matrix bracket is relatively simple, due to the Grassmann charges. 
Consider, for example the case of two particles. The particle or Lienard-Wiechert parts of 
the matrix bracket vanish since Qi = = Q2 ^^'^ cross terms vanish because they involve 
Poisson brackets of particle one variables with particle two variables. Thus the only part of 
the 6x6 matrix bracket that contributes is from the field variables. It has the form 

{Xi(r, ai), X2(r, a^)} = (I - ^)5'(^i - '^2), (5.31) 
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and since 

{Xi,Xi} = = {x2,X2}, (5.32) 

only the 3x3 off diagonal portion contributes. 

In order to have a well defined Dirac bracket we need to use a modified form of the Dirac 
bracket in which the inverse of the matrix of constraint Poisson brackets is used. Calling 
this matrix C, we define so that CC~^ = (I ~ ~ ^2). But the transverse form 

of the delta function allows us to use the idempotent property of the projector to show that 
the inverse of C in this sense is just C itself. In that case for two functions f{Ki, r)i),g{Ki, ffi) 
of the particle variables the Dirac bracket becomes 

{f,9r = {f.9}- 

■' i j 

-{/> -Y^Qj^±Sj{ff - rfji^), ■ {-^QiA±si{^ - Vi{r), i^i{r)),g}]. (5.33) 

j i 

This bracket will lead to a new symplectic manifold by altering the basic commutation 
relations and providing us with new canonical variables. Toward this end we define the 
following scalar function. 

N-l N 

}C=J2 ^ QiQj^iji'^i^i^v'ni-'nj)^ (5.34) 

i=l j=i+l 

in which 

Kij = j £a[A±si{^ - Vi, Ki) ■ 7r±sj{a - ffj, kj) - 

— Hj'i f]i ^j) — A^ji. (5.35) 

Let ffi = ffi + di, Ri = Ki + Pi, i = 1, 2, ..N where 

AT i-l 

di^ tti ^ QiQi^Ki^ij + ^ QiQj^ Ki^ji, 

j=i+l j=l 

N _ i-l 

Pi^bi + h E Q.Q, V„/C,,. (5.36) 

3=i+l 3=1 
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Since they do not appear in these equations, we may choose ai — bi — qn — bN — 0. 

We determine relations between the unknown coefficients by requiring that ffi, Kj be 
independent canonical variables. So, for example, (for k <l) 

{'nk,r}i}* = {Vk,r}i} - 

- [ j d^a{f)k, - Yl Qi^±Si{^ - ffi, Ki)} ■ {- J2 Qj^±Sj{^ - ffj, Kj):fli} - 

i j 

-{vk, - Yl Qj^±sj{^ - ffj, Kj)} ■ {- Y Qi^±si{^ - ffi, Ki),ffi}] = 

= {ffk, + {ttfc, ffi] + QkQN KiK^ki = 0. (5.37) 
Then using the expressions for ctj leads to 

— Ofe = 1, k > I; ai — Uk — —1, / > k. (5.38) 

Solving this gives 

a2 — as — .. — un '■— a, 

ai — a2 — ■■O'N-i '■— o. (5.39) 
Similarly, requiring that = leads to 

&2 = &3 = ■■ = bN := b, 

6i = 62 = .. = feiv-i := b. (5.40) 

^ 

Requiring that {ffi, Ki}* — 1 leads to 

Qi + bi^O, i ^1,..,N -1; CLi + bi^O, i^2,..,N. (5.41) 

This condition implies that the first two conditions are equivalent to one another. The 
requirement that {ffk, ki}* = for A; < / leads to 

ak + bi = l; k<l, (5.42) 



or 



55 



Ofc — a« = 1; k < I, (5.43) 

which is also the same as the first condition. While for / < it leads to 

ak + bi = -1; / < k, (5.44) 

or 

o^k — CLi = —1; / < k, (5.45) 

which again is the same as the first condition. This leaves us with just two unknowns a 
and a. 

So in summary we have 

ffi = ffi + a QiQjVf^Xij + dY^ QiQjV^Xji, 

j=i+l j=l 

Ki = Ki- a ^ QiQjVri.JCij - a^QiQjVri.^ji. (5.46) 

j=i+l 3=1 

Let us rewrite the rest frame condition Eq. (|3.4CI|) 

Hp = V{int) =Y.^i + J d^(^[^± X B]{t, a) = 

i—l 

N 

= ^ + ^ QiQj / d^a Tf^siid" - ffi, Ki) X (V<x x A±sjid' - Vj, %)) + 

j — 1 i<j 

+T^±Sji^ - Vj, f^j) X (Va X A_L5i(a - ffi, Ki))] ^ 0. (5.47) 

in these new canonical variables. Let us expand the cross products, integrate by parts and 
use the transverse gauge condition to get 

N 

rtp = V{int) = E '^^ + E Q^Q3^V,^^J = 0- (5.48) 

If we choose a = —d =\ this becomes 

Hp = V{int) = ^+ = E = 0' (5-49) 

i=l 

like in the case of either free or interacting particles on the Wigner hyperplane [in an 
instant form of dynamics the Poincare generators V^int) do not depend on the interaction]. 
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In other words the rest frame condition is simply that the sum of the N (new) canonical 
momentum is zero. Note that this same choice gives 

AT _ AT 

"^Vi = "^^1 -Yl QiQj'^fij^ij- (5-50) 

i=l i—1 i<j 

Therefore, with the choice a = —a = |, Eq.( |5.46| ) defining the final canonical variables 
becomes 

Vi = Vi+ 2^ QiQi^ i^j^ij^ 

l^i = l^i 2 QiQi^rii^iji (5.51) 

with ICij given by Eq.( p.35| ). 

In the next section we will re-express the other internal Poincare generators M = P("j„t), 
J(int)i ^{int) of Eqs.( [4.1|) and the internal center-of-mass coordinate g+ ~ —IC(^int)/M of 
Eq.( [4.12D in these final canonical variables. 
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VI. THE EXACT DARWIN HAMILTONIAN FROM THE INVARIANT MASS. 



In this Section our aim is to use the exphcit semiclassical Lienard Wiechert solution 
of Eqs.( ^.27 ), ( p.28| ) for the transverse vector potential to obtain an explicit form of the 
instantaneous act ion- at- a- distance potentials present in the invariant mass M = 'P(^„() of 



Eq.(|]T) after the elimination of the electromagnetic degrees of freedom. In its phase space 
form this Hamiltonian for the r = Tg-evolution will contain: 

i) a vector potential 1^(t) = Qi J^l^f Qj^±Sj{Vi{T) — Vji^)^ ^ji'^))-: minimally coupled to 
K,i{T), under the square root kinetic energy term of each particle 'i'; 

ii) a scalar potential f/(r) = ^ / d^o-[T^j_s + (r, a), coming from the field energy, which 
adds to the Coulomb potentials. 

Due to = we can extract the vector potentials from the square roots: the 
semiclassical contribution from the vector potentials is a new effective scalar potential 
Ui = —Y^fLi Qi n-^^ and we get the complete Darwin potential Voar = U + Ui added 
to the Coulomb one. In this form the invariant mass becomes the exact semiclassical Dar- 
win Hamiltonian and Vdar is the Darwin potential to all orders of for every A^. If 
we call Vlod the lowest 1/c^ order historical Darwin potential [see Eq. (|6.11|) j, we have 
Ui = 2Vlod + UiHOD [see Eq.(|JOD], U = -Vlod - U^hod + Uhod [see Eq.(|03)], 



Vdar = Vlod + Uhod- [''LOD" and ''HOD" mean lowest and higher order in 1/c^ re- 
spectively]. When we re-express the invariant mass in terms of the final canonical variables, 
there is an extra contribution U'^od coming from the square roots, so that at the end the 
final Darwin potential is 

Vdar = Vdar + U'^qd = ^lod + ynoD] ynoD = Uhod + U'hod- (6-1) 

In this Section we shall evaluate the complete (to all orders in 1/c^) Darwin potential in 
the old (no longer canonical) variables and then we shall re-express it in terms of the new 
final canonical variables. We begin by obtaining the contribution of the field energy integrals, 
expressed in terms of the canonical particle variables. Using the truncation properties of 
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the Grassmann charges we extract from the kinetic piece the vector potential portion and 
combine it with the field energy integral. In addition to the naive kinetic part (expressed 
in terms of the old canonical momentum) we obtain the rest frame Coulomb part, a part 
that generalizes the standard Darwin interaction and a double infinite series containing all 
higher order corrections. As an extra check on these generalized Darwin interactions in 
M we obtain an independent derivation of this series using the Lagrangian expression for 
the invariant mass. Then we express the kinetic portion in terms of the new canonical 
momentum to obtain the final form of the complete Darwin Hamiltonian. This interaction 
Hamiltonian contains no A^-body forces and is a sum of two-body portions. In the center- 
of-mass rest frame the double infinite series can be summed exactly to obtain a closed form 
expression for the special case of two particles. Then all the generators of the "internal" 
Poincaire' algebra and the energy-momentum tensor are expressed in the new variables in the 
A^-particle case. As with the three-momentum we find that the internal angular momentum 
does not depend on the interaction. Also we obtain the "internal" center of mass g+ and 
there are some comments on how to find a collective variable (replacing the center of mass) 
for a cluster of particles interacting with the remaining ones of the isolated system. 



A. Field Energy and Momentum Integrals. 



Although we have summed exactly to a closed form (|5.28|) the semiclassical Lienard- 
Wiechert solution, we use the series form ( ^.2?| ) for finding the expression for the invariant 
mass M, since the closed form provides no simplification in obtaining that expression. From 
the above we can find expressions for the semiclassical Lienard-Wiechert electric and mag- 
netic fields. For the electric field we find [fjiij) does not contribute due to Eq. ([4.23|) , and 



the same is true of Ki{T)] -^\a — i7j(r)|" = —fj^ir) ■ Vo-|(? — ^i{T~) 
E±s{t, <J) = nj_s[r, cr) = — 



dr 

N 

2m-l 



iV oo -| 
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(2m + 2)! 

N 

^QiE±sti^ - Viir), Ki{T)) = 

i=l 

-?==== A^Sii^ - Viir), Ki{T)) 
i=i ^mf + Kf{T) 

N 
i=l 



An\a - rf,iT)\n |a - r/,(r)| [mf + (/5,(r) ■ ^E|g]|)2]3/2 



The magnetic field is [5^ = e''''^^A\s\ 



I..N 



Bs{t, a) = ^ QiBsi{a - f]i{T), ^^(r)) = 

i 

= - E^ E T^^'.M X V.(,^.(r) . V.)-) |a-r/.(r)|-- = 



In these expressions fJiij) may be replaced with K,i{r)/ymf + R^ij). We have also given the 
closed form of the fields. 

Let us remark that the Feynman- Wheeler complete absorber assumption is violated in 
the Maxwell theory, since it would imply (see for instance Ref. [0, where there is the 
definition of radiation in the Feynman- Wheeler theory) £'_|_5(r, a) = £'_|_±(r, a) = and 
Bsij, a) = B±{t, a) = everywhere (inside and outside the absorbers). 

For the energy we need Ej_g + Bg and for the momentum we need E±s x Bs- In Appendix 
B we evaluate these and show that [fJijiT) = rjij{T)fjij{T) = ffii^r) — ffj^r)] Vij = d/dffij] 

i<j 
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U{r) = \jdMEls + Bl){r.a) := E ^1.) 



and 



An 



iij 



(2n + 2m + 4)! 



271 V72 2n^2m+l 



(2n + 2m + 2)! 



(2n + 2m + 2)! 



(6.4) 



L.iV 



/ d^a{E_i_s X Bs){T,a) := ^ QiQjhi{f)i,f)j,f)ij) 



i<i - m=,n=o (2n + 2m + 2)! 



(2n + 2m + 4)! 



]) + (i ^ j) 



(6.5) 



B. The complete Darwin Hamiltonian in terms of the old canonical variables 

Our first aim is to express the invariant mass M, i.e. the Hamiltonian for the r = Ts- 
evolution, in terms of the original canonical variables. Later we will obtain this Hamiltonian 
in terms of the new canonical variables. The Hamiltonian M and the internal 3-momentum 



are 



i=l 

N 1..N ^ ^ 

^ I CO- 

= E V^f + + E , I : ^ ^ , ,1 + yDAR{T\ 
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VoARir) 



= Vlod(t") + Uhod{j), 



Vi^int) = Hpir) = k+{t) + / £a[7r^s x Bs]iT,a) ^ 



(6.6) 



The first line of M is 



N 



N 



i=l 
N 



i=l 



i=l 



Ki{T) ■ Q^A±s{r,ffi{T)) 



mi + Ki[T) 



'^1 = - 2^ 



1 Jm^ + k2 



VDARir) = f/i(r) + U{r) = VLooir) + UnoDir), (6.7) 



in wliich the vector potential is given by the semiclassical Lienard-Wiechert transverse po- 
tential Eq. (|5.27|) . By re-expressing this transverse vector potential in terms of the momenta 
[one may use in that expression either the new or old canonical variables because of the 
Grassmann truncation] we obtain 

N 



iV CO 1 



V.)^"^|a-rf, 



\2m-l 



so that for the scalar potential Ui we get 



1^1 = - 2^ 



1=1 Jmf + i=i 



1..N 



QiQj 



\2m 



m2 + K,? 



v.,)( 



V,; 



■ V, 



(2m)! 



\2»n+l 



+ ( 



+ 



+ K^ 



V.,)( 



y ^,2^n+l'\ % 
2 ^-^ 



2m+l 



(2m + 2)! 



(6.8) 



(6.9) 
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The lowest order part of this kinetic contribution (the m = term) is twice the famihar 



Darwin interaction (but with rrii — > ^Jinf + nf] strictly speaking this is a higher order 
correction). Thus we have 

N 7^. .y. 

U, = -J2^^^ = 2VLODir) + U^HOD, (6.10) 



=1 + z?,- 



with 



^ ■ V,,)( , ■ Vi,)^) 



, ■ , + ( , ■ r/,,)(^^= ■ r/,,)). (6.11) 



The standard form for the historical Darwin term is the above but with ym? + — >• m^. 

The remaining compensating part of the familiar Darwin interaction plus all higher order 
parts come from the field energy f/(r). In terms of momentum variables, using Grassmann 
truncations, the expression Eq.( |6.4| ) for the field energy integral simply ) becomes 

2^ ^<j 47r Jm^ + K^ J^j + 



47r ±r'.±^, U^2 , -2 /^2 , -2 



i<i m=On=0 ymf + Kf ^TTl^ + 

Kj rr N2m+1/ i^j fj \2n+l 2n+2m+l 

(2ra + 2m + 2)! 

V7 '\2m+2/ «j VT \2n+2„2n+2m+3 

(2n + 2m + 4)! 

(' X7 \2m( ^_^3__ VT N2n 2n+2r?i-l 



m2 + ^2 Jm] + /«| (2n + 2m)! 



\2m+l/ ft^j N2n+l„2n+2m+l 



(6.12) 



(2n + 2m + 2)! 

Single infinite sum pieces can be split off from the double infinite sum in the last two 
lines of the above expression for the field energy integral. This naturally separates out the 
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compensating portion of the familiar lowest order Darwin parts plus all remaining higher 
order Darwin parts, including a piece that cancels exactly Uihod- 

— * — * 



Vij) = 



y^ QiQ, 



1 



K,7 



K, 



N 



o o °° °° 

i<j m=0 n=0 



V7 \2n+l„2n+2m+l 



(2n + 2m + 2)! 



( \7 \2m+2( f^j VT \2n+2„2n+2m+3 



(2n + 2m + 4)! 



+ 



^ 47r ^ 



\ (2m)! 



AC,; 



m,^ 4- Kf 



1 / K 



(2m + 2)!V /^2 + ^2 



• V 



+ 



-V-.f _ ■ V..)2m+1\ 2r 



2m+l 



+ 



I..N 



^ QiQj r ' 

i<j m=on=o y'mf + k| ^Jrr^'^^ 



j 



X 



2m+2/ Kj 



vy \2n+2„2n+2m+3 
V 'Hi 



(2n + 2m + 4)! 



y ^ 2m+3 / «j 



Y7 \2n+3^2n+2m+5 
2 '/ii 



(2n + 2m + 6)! 



—VlOD — UiHOD + Uhod, 



VDARir) = U,{t) + U{t) = Vlod{t) + U„od{t) 



(6.13) 



In this form we see that the first summation is —Vlod and combines with the gauge part 
of the kinetic piece [2Vlod{t)] to give the familiar lowest order Darwin piece Vlod{t)- The 
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third set of summations is —U\hod and exactly cancels with the corresponding term in the 
kinetic piece. The second and fourth set of summations ( the double sums which begin at 
higher order in (1/c^)) we define as Uhod, coming only from U{t). Altogether we obtain 

N 



i=l 
1..N 



N \..N ^ 



1=1 

N 1..N 



E H + ^Kr) + E .... ^fl. + V,ouir) + VHonir) 

i=l i^j I /*\ / /j \ / 1 



AT 



y^;^: ^i^^^m, 21 



i<j 47r „=on=o \Jmj + ^Jmf+^ 

( Kj \J \2m.+l( «j VT \2n+l„2n+2m+l 

(2n + 2m + 2)! 

^ K,: V7 \2m+2/ V7 \2n+2 2n+2m+3 



(2n + 2m + 4) 

vy \2m+2/ \2n+2„2n+2m+3 



^mf + /5| V^T^ (2n + 2m + 4)! 

Kj \7 \2m+3( Kj ^ \2n+3„2n+2m+5 



(2n + 2m + 6)! 



(6.14) 



In Appendix C we show how the multiple directional derivatives in the generalized higher 
order Darwin interactions can be evaluated in general, thereby obtaining a more readily us- 
able form. (The expression for M and for Vdar for arbitrary N is an immediate generalization 
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of the case for N = 2 and for simphcity of notation the details in the appendices will be 
limited to the two-body case.) However, because of the significant complexity of the above 
form it will be of value to first obtain an alternative derivation of this series. For N=2 we 
have the following Lagrangian expression for the invariant mass [see Eq.( [4.10D with A(r) = 0; 
hi is defined in Eq.(p7 



z? hf^ - mi ms Q1Q2 , Q1Q2, ,^ A ^^ 

: = ho{f]i,f]2,f]) + ^^hi{f]i,f]2,f]). (6.15) 
An 

In order to find the Hamiltonian H{k,i, k,2]V) from hi we must demand that Hamilton's 
equation be satisfied . We use the Dirac bracket since we have used the constraint as a 
strong condition on the dynamical variables. This will lead to a set of differential equations 
for H{k,i, K2', v) ^'^d to a result that agrees exactly with the above expression for M. The 
details of this analysis are given in Appendix D. 

Here we mention a comparison of these cross checked results (valid to all order of 1/c^) 



with approximate results obtained elsewhere. In [^ one obtains a single time Lagrangian 
by expanding the symmetric Green function in the Fokker action, used by Wheeler and 
Feynman, to all orders in From this Lagrangian one obtains the Legendre Hamiltonian 

h{fii,r]2, f])+ terms involving the acceleration and all higher order derivatives. Ignoring those 
higher order accelerations one finds the same Legendre Hamiltonian as above. From that 
Hamiltonian the authors obtain a final Hamiltonian that although agreeing through order 
1/c^ with the results above (including the standard Darwin interaction to that order) they 
differ from our common results (above and in Appendix D) at order (no terms of 

higher order are computed in |^). The failure to obtain results that agree with our result 
here is the neglect there of using the proper Dirac brackets in the Kerner reduction. The 
reason that those Dirac brackets were not used is that the authors took as a starting point 
the Fokker action, not taking into account that this action itself is a result of imposing 
constraints on the solutions of the electromagnetic field equations. 
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By using Eqs. (p.51|) we get 



\ 



2'dmi + 



with 



-Axsii^i, Kj, ffi - ffj) ■ {Ri ■ Vf^^TT±si{f^h - ^j))]^ (6-16) 



so that M = becomes (due to Grassmann truncation we can replace the old variables 
by the new variables in the interaction terms) 



N 



N 



(int) / ^ y ■ ■ "« ' ■ ' / ^ j 



+ E 



^'^\ffi- ffj 



mj + K- 



rjil 
2 



+ 



l..Ar 



QiQj 



X ^ '^i^j \ ^ \ ^ r '^j 't'j 

A-rr ^ ^ [ / ~ =2 ■ / 2,-2 



2n+l ~2n+2m+l 



(2n + 2m + 2)! 



'\2m+2f' vy \2n+2 ~2n+2m+3 



{2n + 2m + A)\ 



+ 



Y7 \2m+2( «j Y7 \2n+2 ~2n+2m+3 



■2 ^-^ ' 



(2n + 2m + 4)! 



YJ ^2n+3 ~2n+2m+5 



(6.17) 

(2n + 2m + 6)! J ^ ^ 

where = d/dfjij. We emphasize at this point that the interaction terms are all two-body 
forces; there are no body forces in our semiclassical treatment. 
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C. Complete Darwin Hamiltonian in terms of the Final Canonical Variables for the 

N-Body Problem. 



In this subsection we will present the final expression for M = VJ^^^^^ in the final canonical 
variables for arbitrary . Now from Eg. ( ^.48|) above we can replace J2i<j r,iQiQj^ij by 



— / (i'^o-[7r_L5 X Bs\{T,(y). We have already performed that integral (see Eq.( pi4| )) for the 
two-body problem and the result for the body problem is an immediate generalization. 
Using that complete expression we obtain [rjij = \ fiij \ = — r]j\] 



N 



i=l 



Qi^^j V — ^ V — ^ r 1 



2m ( T-T \2n+2\ ~2n+2m+l 



+ 9f — . v:.)2m+2/ _ . y..N2n+2 , 



■ V:,)^'"+^(^i^ ■ Vl.O'-^^'jC'"^^' • (6-18) 



We now combine this portion of the complete Darwin Hamiltonian coming from the vector 
potentials Viir') in the kinetic terms with the potential Vir') arising from the field energy 
integral [which contains VLOH{jy\ to obtain the final form of the invariant mass in the rest 
frame 
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M . 

i=l 



N 



QiQj 



i=i i<j ^^'Hj 



N 



QiQ j 



E V + (^) + E + ^i^AK(r) = [Vdar = Vlod + Vhod] 



i=l i<j iV 



Vdar — — E 



QiQj 



iA 47r 



, ^ ^' / ^ ^-VyA I ^ 2 



QiQj sr-^ sr-^ f l^i K,j 1 



8^ 



~2n+2m+l 



I / '''' xf VT* \2n+3\ --2n+2m+3l 

I ^— :^-'^ij) \ I Z^'^^i) J'lij J 



( iij vy \2m+2/ \2n+2 ~2n+2m+3 

I / _ _2 ■ ^ij^ V / , ~2 ■ ^V) 'lij 



_ l^i l^j V "'H^i 

4^.' (2n + 2r« + 4)! 



V7 N2m+3/^_K2^_ Y7 \2n+3^2n+2m+5 



+ (2n + 2,„ + 6)! \- («^^^* 

Notice that the N body interaction Hamiltonian is a sum of individual 2 body Hamilto- 
nians.(Grassmann truncation eUminates A'"-body forces). In this sense these interactions 
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correspond to a sum of disconnected and spectator Feynman diagrams. 

D. Further Reductions of the N-Body Darwin Hamiltonian and Closed Form 
Solutions for the Two-Body Problem. 



It is of interest to see if we can make further simphfications by deriving an expression 
for the multiple derivatives. We show in Appendix C that the problematic expression (/tj • 
^ijYi^j ' Vij)^f/'^'^^~^ is obtained in terms of powers of [fiij = |r/j — f]j\; fjij = {fji — f]j)/fiij] 

COS^ 0,, := {^^■^f^J-f^^■V^p■ViJf _ ^g_20) 

For N = 2 this expression reduces to 1 if ki = — ^2 (the center of mass rest frame condition 
for the two body problem) . The result for a = 2m + 1, 6 = 2n + 1 is 

[(2m + 2n + l)!!]2[(n + m + 1)!]2 
(2(m + n + 1))! 

(Kj ■ Kj — Ki ■ fjijKj ■ flij){^K.i — {t^i ■ flij)'^)"^{K.j ~ {t^j ■ fjij)'^)^ 

~ X 

Sis (' :'::-«) P" ' &-^rH.^'-^-- 

For the two body case in the center-of-mass rest frame ki = —k.2 '■= this expression 
reduces to 

- l(!!!i±i!i±i)!!!!,g^ _ (g . (6.22) 

V 

In the general case this expression can be more simply written in terms of cos(2A; + l)4>ij 
[see Eq. (|C55|) ]. Using the identity 

[(2m + 2n + l)!!][(n + m + 1)!] 1 



(2(m + n + l))! 2"^+"+i^ 
one obtains for m > n the expression 



(6.23) 
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2 (2(m + n + - {k, ■ f^^,fr+^'\K] - • 

(2n + l\ /2m + 1 



X 



Ecos(2A; + l)0,,( J( J (6.24) 



A;=0 

Hence 



\n — k ) \m — k 



2m+3/-t V7 N2n+l~2(m+n)+l 



2 (2(m + n + 2))\{k' - {k, ■ - {k, ■ %)=^)"+V^ 



fe=0 

By similar methods one finds that 



fjij 2"^+"+4 

/x^ .X, /2n + l\/2m + l\ /2m + 2\ /2n + 2\ 

To determine the other combination note 



2 m:n< + n))K^-(i:!rr,,j)r'(^-(^j-nu)r ^, 

fj.. 2^+'*+2 

/'^^ / , X , / 2n - 1 \ / 2m - 1 \ /2m\ f2n\ 

cos(2* + 1)*, L _ 1 _ J L _ 1 _ J + L L 1 (6-27) 



fe=0 

Thus 
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2 (2(m + n + 1))\{k] - ■ ^l^,fT+\K^ - (k, ■ ff^^fY 
fjij 2™+"+3 

/ , , , / 2n - 1 \ / 2m + 1 \ (2m + 2\ (2n\ 

Unfortunately, the A; sums cannot be performed analytically to known closed forms for the 
A^-body problem. 

In the special case of the two body system we can obtain a closed form if we use the 
rest frame condition /?i + /?2 = . The expression we get in this way may be used with 
the Dirac brackets associated with ki + /t2 ~ 0, ~ 0, so that the final reduced phase 
contains only fj = — 172I and /?:=/?! = —hi2- Using the identity in Eq. (|6.22|) the higher 
order Darwin part Vhod = Uhod + U'^od becomes 



(2n + 2m + 2)! ^ ^ ^ 

x( , ^ , ^ + — ^) + 



[(2m + 2n + 3)!!]2 



(2n + 2m + 4)! 

t2 I ,-f ^ hm'i I m.-r -I- ^r, rr7,„ 



'mf + R \lml + K mf + K + k 

^2 [(2m + 2n + 3)!!p ~2 _ ^ _ ^l2)n+n^+2( 1 )2n^+3( 1 )2n+3 



mf + K^ \Jml + K 

[{2m + 2n + 5)!!]^ ^^2 ~ ^.^2in+m+3^ ^ ^2m+3^ ^ ^2n+31 

(2n + 2m + 6)! " ' ^ ^X^^^ ^XT^' 



mf + K Y m2 + K 



(6.29) 



We use 

[(2m + 2n + 1)!!]2 _ (_)«+™ /-3/2' 



(6.30) 



(2n + 2m + 2)! 2{n + m + l)\n + m) 
and let m + n = / so that < m < / and < / < 00. Then we perform the m sum using 



m=0 

and obtain 



1// y{y-x) 
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1 ^21+2 



Vhod — 



Qlg2y^^_^2 (- 



2(/ + l)V / 



•^21+2 _ ^ 



J ■\2l+2 



mr — m5 



^ + ^) + 

2 1-* 2 1-* 

mf + K, 777.2 + 



J \2l+2 



2{l + 2)\l + l) 



yi+2 _ ^ 



1 \2l+2 



9 



2 2 
777f — 7772 



xy7nf + K \lm2 + K 



~2 ' ~2 
7?7f + K, 7772 + K 



^2) + 



J \2W 



-)2/+2 _ ^ 



nn\+K 



2 2 
777i — 7772 



-]( 



1 



r)(- 



1 



)- 



(/ + 3) + 2 

\jm\+i' \Jml+ft 

ml — ml 



r)(- 



(6.32) 



We use 



E 



^2(; + i)V / 

00 (' V+i 



E 



S 2(^ + 2) 

^ V~/ /^~3/2^^;_ 1 



a;)-3), 



a; 



^2(/ + 3)l,i + 2j 
and finally determine 

Q1Q2 



15. 
8 ' 



(6.33) 
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ml + /S? 
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+ ~2 '^J 



-{ml + K)[2^ 
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2 , -2 

777f + «; 



+ 



\^ml + {K- 77)2 ^2 ^ 
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X 
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2 , -2 

mi + K 
— 5 + 



\ ml + K \ ml + K 



7712 + 



) + 



>2 + '? )[2, 



~2 ~ 
+ K 771^ + (^-77)^ 

+ ~2 "^J 



^ 777i + • 7^)2 ^2 ^ - 



-(777? + [2 



77lf + K 



^)( 



\Jml + K 



\^ml + {K-f]f ml + K 

' )- 



X 



7772 + V 



(7772 + )^ 



777; 



+ 3 / 7772 + ('^ ■ VY 



\ 777i + • 77)2 8 V mi + 



+ 



5 1 ml + (/v • iif 

^ \ 7772 + K 
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+ 



+ 



lml + (K- 7/)2\ 




\ 777f + K / 


8 



X 



(- 



f)(- 



\Jm\-\-K' \jm\^K 
So our final two-body expression is 



M 



777? + ^' + V777i + ^'+^ + W 

Vdar = ^ozi + Vhod, Vhod = Uhod + U'jjod 



(6.34) 



r)( 



'7r7f + \/ 7772 + 

in which we have used the rest-frame condition ki — —K2 '■— Note that Q1Q2K — QiQ2K- 
In the equal mass limit 7771 —> 7772 := 777 



)■ 



(6.35) 



VHOD — H — ^—r- X 



m^[3ti + ■ fif] - 2k^[k -3(/?-7))2] 



m2+(K-r))2 



[3^t (k • 7))2][7772 -f- {K-fjf] 



(7772 -|- «;2)[j7^2 _|_ ^j^ . ^^2j 
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(6.36) 

so that 

i2 , Q1Q2 . Q1Q2 



M = V;,^,^ = 2Vm2 + K + ^ + ^ X 



m^lsi + {k ■ fff] - 2i[i - 3{k ■ fj)\l f+f -2[i + {K- r/)2][m2 + (/? ■ r))^] 

(m^ + K ) [m^ + (k ■ r))^] 

(6.37) 



Our result differs from tfie appropriate terms at order with those in |^ and ||50 
although it does agree with their Darwin interaction at order Since each of these 

latter three sources took as a starting point the Fokker particle Lagrangian they have not 
used, as we have done here, the canonical variables which would come from using the pair of 
secondary constraints arising from the reduction of the field plus particle Lagrangian to the 



Fokker action. We point out that Molina et al in [^, like earlier work by Golubenkov and 



Smorodinski |^obtain order 1/c^ corrections to the Hamiltonian. However, unlike them 



these authors not only use the Coulomb force law to replace acceleration dependent terms 
on the Lagrangian but also include the effects of that substitution on the choice of canonical 
variables by viewing it as a constraint (thus including Dirac brackets), in the reduction 
to Hamiltonian forms. (Note that since these three approaches, unlike ours, do not use 
Grassmann charges, they contain acceleration driven terms not contained in our approach. 
The comparisons we are talking about here with our approach refer to the terms not driven 
by acceleration dependent Lagrangian potentials). 

It is of interest that Hamilton's equations has a solution for circular orbits just as does 



the Schild solution for Feynman- Wheeler electrodynamics |^ . In Appendix E we show how 



this comes about for equal masses. The case of unequal masses is similar. 
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^* '^{int)^ ^\int)i ^+ Energy-Momentum Tensor in the Final Canonical 



Variables. 



Eqs.(|^ and (Eq.(|6l9D [(Eq. (|05|) ) for = 2] give V(^int) and M = P(^.„j) in terms of 
the final canonical variables. For the internal angular momentum we get 

i=i •' 

i=l 

1 ^ ^ r - - 



(6.38) 



Using the forms for cti and (3i given in Eqs.( |5.5lD together with the expression for ICij and 
expanding the cross products in the integral we obtain 



N 



JLt) = Y.ivi X KiY + 



1=1 



1 ^ ^ /■ - 

+ 9 QiQjiiVi X Vr,, + X V^ J / d^a{A±si ■ T^LSJ - A±sj ■ T^±Si)} 

1 ^ ^ r 



1=1 j^i 



(6.39) 



Using the transverse nature of the field together with vanishing surface terms we find that 



N 



J{int) =^rji^f^i 



(6.40) 



if 



{Vi X V^, + X VkJ j d a(A_L5i ■ vr^sj - ^±5i ■ Tr^^i) 

d^(T[A^St X 7fx5i - X (V^AV5i)7risj- + A^sj x vf^s* + ^ x (AVsiV^vrl^^)]. (6.41) 



We have given explicit forms forms for A^si and T^±si in Eqs.( ^.28D and (|6.2| ) respectively. 
Their general forms are 
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^±Si = -p-R- — ^[i^ifiii^'i, i^i ■ pi) + Pi9i{i4, i^i ■ Pi)], 
A7T\a - r]i\ 

T^±Si = ■ 1^"^ ^ [Kihi{Kj, Ki ■ pi) + piCi{Kl, Ki ■ Pi)]. (6.42) 

in which Pi := -fe^ . Usiiiff 

\(J-rii\ 

V.^^ = V.p. = -Zl^(I-p,p,), (6.43) 

we find that 

[{m X + X V^Jv4x5i] ■ 7rx5i = -<?>< (V^v4^±5i)7rl5j- + A^si x ^±5^, (6.44) 

while 

- A^^SjM X V^, + /?i X V^Jvri^^ = A^xsj*? X V^-n'^^si + ^±5i x vf^^i, (6.45) 

thus verifying Eg. ( |6.40|) . 

Thus as expected, in an instant form of dynamics, the internal angular momentum does 
not depend upon the interaction. On the other hand, the interaction-dependent internal 
boosts have the form of Eq.( [4.1|) . Using the above forms for the final dynamical variables 
we obtain 



N 



^{int) = - H 'Ti [ V ITT'i +f^i + 
i=l 

f^i ■ Y.jjtiQiQj[^fi,lK,ij{f^i, f^j,Vi - ffj) - 2A±sjifij,'ni - ffj)]^ 



~9 QiQj\l^l + /«/VsXij(Ki, Kj,ffi - f]j) + 

QiQi ffj — ffj ^ y, „ QiQj f ^s^T^lsji^ " Vjj i^j) 



1=1 j^i 

N 



\Vi-Vj\ i=ij^i 4vr J \a-f]: 



~2 ^ ^ ^'^^ J ^^^^I'^^Sii^ " ^i) ' ^-L5j('? - Vj^ + 
i=l j^i 

+Bsi{a - ffi, Ri) ■ Bsj{a - ffj, fij)] = 

-Vlir,t)R^. (6.46) 
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In the last line the internal M0ller center of energy ( [4.12|) is shown explicitly. This equation 



allows us to express the internal canonical center of mass gV ~ i?+ defined in Eq.( [4.12D in 
terms of the final canonical variables. The natural gauge fixing to the rest frame constraints 
V{int) = Eili ^ is g+ ^ 0. 

From Eq.( pj.44D we get the following expression for the conserved energy-momentum 
tensor: 



N 



i=l V jV« 



1 d d 



Bsi ■ Bsj]{Ts, a), 



A 



i=l jyti 



1..N 



d 



N 



T-(T„a) = E'5'('?-^.m)) 

i=l 



1 N -I f) f) 



+ Bs, ■ Bsj] 



QV 



gs 



+ B^SrBU]iT,,a). 



(6.47) 



By using Eqs.( |5.51| ) we get [r/^ = r/j(T^), Ki = Ki{Ts); A±si{a-f]i, Ki) is given in Eq.(^ 
iCij = IC,j{Ki, Kj,f]i - ffj) is given in Eqs. (|5.34|) , ( |5.35|) ] 



); 



TV 



1..N 



T^^{T,, a) = E s'i^ - + ^ E Q.Q. v.,/c.,) 
i=l ^ j^i 



\ 



1..N 



+ 
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1 ^-^ ^ ^ Q 



d 



+ Bsi{a - ffi, Ki) ■ Bsj{a - ff., Kj) 



N 



1=1 



+ 



1 d 



■{E±sj{B - ffj, Kj) + -^5^ (a - f]j) + 
+ Bsi{a - ffi, Ki) ■ Bsj{a - ffj, Kj)] , 



N 



1 



l..Ar 



1=1 j^i 



~< -QiY. Qj {^IsjiVi - Vj, i^j) + ^V^Xij 



+ 



1..N 



xBsj{a -f]j,Kj) + {i ^ j) = 



i=l 



2 

+ 



1..N ^ ^ ^ Q 

+ QiYl l(E±si{a - fji, Ki) + —S\a - fji)) X 
xBsj{a -f]^,Kj) + {i ^ j) 
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N 



1..N 



i=l j^i 



~ ~ ~ \ _ 



1..N 



d 
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■[E^sA^-fj,,K,) + -5\a-fj,)] + 



+ Bsi{B - ffi, Ki) ■ Bsjia - f]j, Kj)) - 
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AT ^ 



1..N 



2 , -2 

1 



Qi E Qj f^li^lsjivi - Vj: ^j) + o v^Xi,) + 



+ 



\..N 



Qi E Qj [b^'([E±sr{B - ~Vi, ~^i) + - ~fi,)] 



ji^i 



d 



+ Bsi{a - fji, ki) ■ Bsj{ff - ffj, kj)) - 

d 



- [Elsii^ - Vi, Ki) + ^5\a - vMElsii^ - ffj: i^j) + ^S'{a - ffj)] 
Esii^ - ffi, fii)Bsj{B - ffj, kj) 



(6.48) 



In the final canonical variables there is a dipole term [gradient of delta function] for 
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each particle like it happens with spinning particles 0: the role of the spin is taken by 

Let us remark that following we can define Dixon's multipolar expansion of 
the energy momentum tensor about an arbitrary point on the Wigner hyperplane. There, 
the requirement of having the mass dipole vanishing could be shown to be equivalent to the 
identification of the point with i?+ namely the internal center of mass. 

For a cluster of n particles inside the isolated A^— body system, we can now define a 
non-conserved energy- momentum tensor T^^{Ts,(j) by collecting in the previous equation 
all the terms depending on the canonical coordinates fji, Ri of the n particles of the cluster. 
This cluster energy-momentum tensor depends also on the canonical coordinates of the other 
N — n particles: for the cluster these are "external fields" . If we make a Dixon multipole 
expansion of this cluster energy-momentum tensor with respect to an arbitrary point on the 
Wigner hyperplane and we require that the cluster mass dipole vanish, then we can identify a 
M0ller center of energy i?!^'* for the cluster. This is the only collective configuration variable 
which can be defined for a non-isolated cluster, which has no internal conserved Poincare 
algebra associated with it, besides the nonconserved cluster 3-momentum V^^^^t) ~ ^ie{n} ^i- 
Moreover, while T^^^{Ts,(t) is the (non-conserved) energy density of the cluster, by analogy 
with the theory of dissipative fiuids we can say that: i) q^iTg, a) = ei^ {u{ps))T^n)i'^s, ^) is 
the heat fiow; ii) V{T„ a) = | Er T[^)iTs, ^) is the pressure; iii) Tl^Yan)iTs, a) = T[^)iTs, a) - 
I J2u ^(n^l^s) ^) is the shear (or anisotropic) stress tensor. 
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VII. CONCLUSIONS. 



In this paper we analyzed how to extract the act ion- at- a- distance interparticle potential 
hidden in the semiclassical Lienard-Wiechert solution of the electromagnetic field equations, 
a subset of the solutions of the equations of motion for the isolated system formed by scalar 
charged particles plus the electromagnetic field. The problem is formulated in the Wigner- 
covariant rest-frame instant form of dynamics, which is defined on the Wigner hyperplanes 
orthognal to the total time-like four-momentum of the isolated system and which requires 
the choice of the sign of the energy of the particles (in this paper we considered only positive 
energies) . 

This was possible due to the semiclassical approximation of using Grassmann-valued 
electric charges (Qf = 0, QiQj 7^ for i 7^ j) as an alternative to the extended electron 
models used for the regularization of the Coulomb self energies. How this happens was 
shown in Ref. |^, where the Coulomb potential was extracted from the electromagnetic 
potential by making the canonical reduction of the electromagnetic gauge freedom via the 
Shanmugadhasan canonical transformation. This is equivalent to the use of a Wigner- 
covariant radiation (or Coulomb) gauge in the rest-frame instant form. 

Ref. [0] presented the retarded Lienard-Wiechert solution for the transverse electromag- 



netic field in the rest frame instant radiation gauge: in this gauge, due to the transversality, 
the retarded Lienard-Wiechert potential associated with each charged particle depends on 
the whole past history of the other particles. At the semiclassical level a single charged 
particle with Grassmann-valued electric charge does not radiate even if it has a non-trivial 
Lienard-Wiechert potential, avoiding therefore the acausal features of the Abraham-Lorentz- 
Dirac equations, and has no mass renormalization. However, a system of A^ charged particles 
produces, by virtue of the interference terms from the various retarded Lienard-Wiechert 
potentials of the particles, a radiation which reproduces the standard Larmor expression 
for radiation in the wave zone, when the particles are considered as external sources of the 
electromagnetic field and their equations of motion are not used. 
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If instead the particles are considered dynamical, the use of their equations of motion and 
of the semiclassical approximation lead to a drastic simplification of the Lienard-Wiechert 
potentials and fileds. Indeed, if we make an equal time expansion of the delay by expressing 
these potentials and fields in terms of particle coordinates, velocities and accelerations of 
every order, it turns out that all the accelerations decouple at the semiclassical level due to 
the particle equations of motion. 

Therefore, at the semiclassical level the retarded, advanced and symmetric Lienard- 
Wiechert potentials and the electric and magnetic fields coincide and depend only on the 
positions and velocities of the particles, so that we can find their phase space expression in 
terms of particle positions and momenta. 

In this way the semiclassical Lienard-Wiechert potential and fields can be reinterpreted 
as scalar and vector interparticle instantaneous action-at-a-distance potentials. It is then 
possible to identify a semiclassical reduced phase space containing only particles by elimi- 
nating the electromagnetic field by adding by hand second class contraints which force the 
transverse potential and electric field canonical variables to coincide with the semiclassical 

— * — * 

Lienard-Wiechert ones in the absence of incoming radiation: A±{T,a) — A±mr{T, a) ~ 0, 
7fx(r, a) — 'ff±LwiT, a) ^ 0. Let us remark that this could be done also in presence of an 
arbitrary incoming radiation A±(^rad){'T, T^±{rad){'T, ^) = — ^^±(ra<i)('?") ^) [it is an arbitrary 
solution of the homogeneous wave equation and must not be interpreted as a pair of canonical 
variables] by modifying the constraints to the form A^{t, a) — ^±Lvr(T, a) — ^_L(rad) i'^-, ^) ~ 0, 
7r±(r, a) - 7r±Lw{r, a) - n±^rad){r, a) pa 0. 

The reduced phase space is obtained by means of the introduction of the Dirac brackets 
associated with these second class constraints. Since the old particle positions and momenta 
are no longer canonical in this reduced phase space, we had to find the new (Darboux) basis 
of particle canonical variables. The generators of the "internal" Poincaire' group inside the 
Wigner hyperplanes in the rest-frame instant form of dynamics can be reexpressed in terms 
of these new variables: the 3-momentum V(int) and the angular momentum J7(int) become 
equal to those for free scalar particles (as expected in an instant form). The interaction 
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dependent boosts /C(j„f) are proportional to the "internal" canonical center of mass g+ inside 
the Wigner hyperplane: gV ~ are the gauge-fixings to be be added to the rest-frame 
conditions V{.int) ~ 0, if one wishes to re-express the dynamics only in terms of particle 
"internal" relative variables. Also the energy-momentum tensor has been evalulated in the 
new canonical variables and there is a suggestion on how to find the M0ller center of energy 
of a cluster of n particles contained in the N particle isolated system. 

The Hamiltonian in the rest frame frame instant form, generating the evolution in the 
rest-frame time of the decoupled "external" canonical center of mass, is the "internal" en- 
ergy generator M = (the invariant mass of the isolated particle system). The 
semiclassical Lienard-Wiechert solution implies the existence of interparticle action-at-a- 
distance potentials of two types: vector potentials minimally coupled to the Wigner spin 1 
particle three-momentum under the square root associated with the kinetic energies; ii) a 
scalar potential (including the Coulomb potential) outside the square roots. In the semi- 
classical approximation all these potentials can be replaced by a unique scalar potential, 
which is the sum of the Coulomb potential and of a generalized Darwin one for arbitrary 
A^. It is the (semiclassical) static and non-static complete potential corresponding to the 
one photon exchange tree Feynman diagrams of scalar electrodynamics and is a completely 
new result. The expression we find contains no A^-body forces, being simply a sum of two 
particle interactions. This is a consequence of our use of Grassmann charges. 

In the N = 2 case we obtain a closed form of the solution by evaluating it in the rest 
frame after the gauge fixing ^ 0: the lowest order in 1/c^ contribution of the generalized 
Darwin potential agrees with the expression of the standard Darwin potential. We then 
show that in a semiclassical sense a special solution of the Hamilton equations is the Schild 
solution in which the two particles move in concentric circular orbits. We evaluate the 
frequency for equal masses 

Future work will proceed along three parallel courses. 

The first is the extention of this work to include semiclassical spinning particles. This 
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will not only build upon the work here but also that of Ref . . Of particular concern there 
will be the issue of reproducing the correct spin-orbit and Darwin terms of the appropriate 
order (the two-body extentions of such terms in the one-body Dirac equation). 

The second line is that of the quantization of our general Hamiltonian. Let us remark 
that quantization of the closed form Hamiltonian for = 2 in configuration space would 
involve not only the usual nonlocal operators for the kinetic energy but also non-local oper- 
ators for the Darwin portions of the potential. For free scalar particles the positive energy 
wave equation idr(p{T, a) = \/m'^ + A0(r, a) has been studied in Ref. Il0| by using pseudo- 



differential operators. Instead see Refs. pj] for the difficulties in quantizing the equal mass 



two-body problem {H = 2\Jvn? + + ^) with only the Coulomb potential outside the 
square root. Note that the so-called spinless Salpeter equation would correspond to the 
quantization of our Hamiltonian with just the Coulomb interaction and at most the lowest 
order Darwin interaction. For the quantization when there are only scalar potentials inside 



the square root (ym^ -|- V^(|f/|) -|- k^) see Ref. p7[ . 

In our semiclassical approximation we have two options for the action-at-a-distance po- 
tentials: i) the vector ones under the square roots and a scalar one V outside; ii) a unique 
effective scalar potential sum of the Coulomb and Darwin (Vda/j = Vlod + Vhod) ones 
outside the square roots. It can be expected that the results of the quantization of the two 
options would produce inequivalent theories. 

The third line is the development of the quantization of scalar electrodynamics on the 
Wigner hyperplanes in the rest-frame instant form. This would be a special Wigner-covariant 
instance of Tomonaga-Schwinger quantum field theory with a well defined covariant concept 
of "equal times". To introduce a particle concept in such a quantum formulation, one 
will have to define Tomonaga-Schwinger asymptotic states and a reduction formalism. The 
natural candidates for the A^-particle wave functions in such asymptotic states in the case of 
the Klein-Gordon field would be the wave functions corresponding to the quantization of 
positive energy scalar particles on the Wigner hyperplanes in the rest-frame instant form. 

Moreover, in the rest-frame quantum field theory there will be new covariant "equal time" 
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Green functions. This should allow the definition of a relativistic Schrodinger equation for 
bound states (replacing the Bethe-Salpeter equation for N=2 and avoiding by construction 
its problems with the spurious solutions, which are a byproduct of the use of asymptotic 
Fock states, since in a tensor product one cannot eliminate the possibility that an "in" 
particle be in the absolute future of another one) . In the rest-frame quantum field theory it 
should also be possible to include bound states among the asymptotic Tomonaga-Schwinger 
states: they should be described by the quantization of isolated N particle systems hke the 
one studied in this paper. 

ACKNOWLEDGMENT: L.Lusanna wishes to thank Prof. K.Kuchar for helpful dis- 
cussions about the need of compatibility with the equations of motion of the second class 
constraints to be added by hand to select a symplectic subspace of solutions in the space of 
solutions of a given isolated system. 
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APPENDIX A: EXACT SUMMATION OF THE TRANSVERSE VECTOR 

POTENTIAL SERIES 



In this Appendix we perform explicitly the summation of the vector potential below [see 



Eq.( |5.27|| with ff^ = Pi replaced by /tj/ym? + Kf] 



2^) 1^ _ ^.|2m-l 

,=1 47r ^0 L (2m) ! + ^2 ^ ^^2 + ^2 ' 

^ / \2m+l I-? --»|2m+l 



: = ^±i(a,r) + ix2(a,r). (Al) 

Using the result of Appendix C that (k^ ■ V^)2"|7/|2'"-i = [(277i - - (k^ ■ 

we get [V^ = d/df]\ 

A ( ^ r [(2777- 1)!!]^ 2 ^-Vi ^2^n^1 ..9^ 

By using 

[(2777- 1)!!]2 (2777)! (777 - 1/2)! 0F(-)™-l (777 - 1/2) 



we find that 



(2777)! (777!)222™ 0F777! (1/2 - 777) !r77! 

(-1/2 -77l)!r77! ^ ^ V 777 y 



(A3) 



i=i^'^ ^ rrii + Ki \(j — rii\ m=o 

= y^, ^ , , , , ■ (A4) 

^^4vr|a-7^,|^^2 + (^^.^p 

For Ax2(t, 0^) we need an expression for (/tj ■ V,,)^™'^^|7fp'^+-'^. One can show by an 
induction procedure that 

m / \ (f^. . ^-Vi ')2;+l 

(/^. • v.)^-+v - mr"-' = [{2m + mr e(-)i 7 ^^'y ^ 21+1 ' ^^^^ 



and hence 
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i=i 4vr ^/ttjS _^ ^. 



E 



[(2m + l)!!]2 



z=o 



(2m + 2)! V"^i+'?i, 



X 



2/ + 1 



(A6) 



Now 



[(2m + 1)!!] 
(2m + 2)! 



^m+l j 

\m + L 



(A7) 



so that 



Ai_2{t, a) 



E 



1 ^TT ./TTi^ + m=0 



E 



-1/-V2' 



i=l 471" ^/ 77^2 _^ ^2 



\ |(T-J7j| ^ 



1 



w 



2\m 



V '^i / 



-E 



Qi V, 



1 . /„i2 + ^2 Jo 



, m + 1 / V mf + K,f 



X 



^m + 1 y y mf + 



) ( — w"^ 



mf + nf 



E( 

m=0 



\r?i+l 



-1/2 \ (nf-w 



2 \ 



E— V, 



,m + 1 / V m,- + 



i=l 



mf + K,f ^ 



1 



+ w"^ 



X 



mf + 



dw. (A8) 



Now 



- r/j 
1^ - fji 



\o- - r]i 



■I 



|a - f]i\ \a - ffi 



■)■ 



(A9) 



So we obtain 



A±2{t, a) 



N 

E 



Q^ 



^ 47r |a - ffi 



mf + (k. 



^12 

* " \^-m\> 



a -rji a-rii 
\a - ffi\ \a - ffi 

1) 



X 



mf + Rf 



(AlO) 



which when combined with the expression for A±i{t, a) yields the results given by Eq.( ^.28| ). 
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APPENDIX B: COMPUTATION OF FIELD ENERGY AND MOMENTUM 

INTEGRALS 



Here we carry out the details in the computation of the field energy and momentum for 
the case N = 2. The general N results obtained in the text are an immediate generalization. 
From Eg ■ ( |6l2D and Eg. ( |6.3|) we find that 



^2 / \ Q1Q2 1 1, 1, 



Q1Q2 



00 00 



Svr^ (2n + 2)!(2m)! 

Q1Q2 1 
"^io^o (2m + 2)!(2n)! 

[iff, ■ V.r^')\a - V2\'''-'][iff, ■ ■ V.)^"^+^)|a - rf,r^'] + 

Q1Q2 1 



m=0 



„t^o(2m + 2)!(2n + 2)! 



Q1Q2 ^ ^ ^ 



^ont^o(2^)K2m)! 

■ V.)(^^ ■ V.)^")|a - r/ir™-^] ■ [{ff, . ■ V.)^")|a - ff,\'-\ (B2) 

Q1Q2 1 u u 

Q1Q2 1 u 

2^ 2^ fo^\uo^\} ^2 X 



16vr2 ^on=o(2n)!(2m)! 
Q1Q2 ^ ^ 



167r2 ^„^j2n)!(2m + 2)! 
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^ Q1Q2 ^ ^ 



\2n-l 



16^' ^on=o(2^)K2m + 2)! 
\a-ffi 

+ — ^2). 



V2 X 



Our aim here is to compute 
l/dV(ii, + 5l)(r,a): ^^^^ 



EE 



1 



167r2 ^o;S.(2^)K2m)! 

"-^2mn 



(^1 • V2)hmn - 
1 



(2n + l)(2n + 2)" 



'3mn 



+ 



1 



(2n + l)(2n + 2)(2m + l)(2m + 2) 



(2m + l)(2m + 2)" 



ldME±sxBs){r,a): 



QiQ 



00 00 



1 



167r2 ^or1^o(2^)K2m)! 
1 ^ 



(2m+ l)(2m + 2) 
+(1-2), 

which involve ten different integrals defined by 



'9mn 



(^1 ■ V2-^7mn V2-^8mn 

V2 



(2m+ l)(2m + 2) 



' lOmn 



— * 

^Smn 
— * 



= / rfV 





— * 

a 






a 






— * 

0" 





V2 



\2n-l 



d^a 

d'a[V^{ff,-V,f^\a-ff2\''''-^ 
Jd^V,{ff,-V,)'^+'\a 



{ff2-^am-^.rn^-m\'"'-' 



(^''2-V.)(^i-V.)^ 

(^\-v.)(^2-v.n^-^2r"- 



^2m+2|^ _ ^^j2m+l 



^2n+2| -; |2n+l 

\cr — r]2\ 



V.(%-V,)^"|a-r/2 



\2n-l 



{Vl-^.){V2-^a)n^-V2 



\2n-l 



iV2-^a){fjl-^.) 



2"i+l|^ _ ^_^|2m-l 



^1 



2m+l 



V.(r;2-V,)^"|a-r;2 



|2n-l 
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The powers of a in each of the ten integrands is cr^"* ^ ^cr^" ^ i ^ ^ 4 Therefore the 
integrals converge. Thus 



lim 



Jo 



da / dila( 



lim 



lim 



(B7) 



A— >ooio J A^ooA A^oo 

Since /j(A) is finite for finite A, we can bring out the derivatives. Now perform integrations 
by parts, change to d/dr] from d/da, translate, and then use the fact that the integrals are 
finite so that they can be replaced by Jj^. This gives (77 := ffu — rfi — 772) 



-^2mn 
Izmn — 
I Aran 
Ibmn — 

— * 
— * 

Iiomn — 



J A 



N2n+2 / 73 



dV|a-(rfi-7f2)|^™+V 



-^'M ■ ^vf'^{f}2 ■ V,f-jfa\a - (771 - 772)1^"^- V^--\ 



\2n+l 



2m-12n-l 



d'a\ff-{7fi-ff2)r-^a 



-^M ■ ^vf"'^\n2 ■ V,)^'^ / d'a\a - (771 - 772)r"^-V2'^-\ 

J A 



im ■ ^,f "'^\v2 ■ V,)^" / d'a\a - (r/i - ff2)r-'a 



2m-12n-l 



J A 



|2m+1^2n-l 



2m+1^2n-l 



The integrals that remain to be evaluated are each of the form 

Ja Jo J-1 

7]{2m +1) Jo ^ 



2^^2n+2m+l 2m+l 2m + 1 

^ 2n + k + l 



(B8) 



■ - 1] - 2 ^ (B9) 

^ k ^ ^ 2n + k + V k ^ ^ 



(2m + 1) 

It is of interest to show how the A dependent terms vanish after the 77 derivatives have acted. 
First note that only the even k terms of the A dependent sum survive. Of the ones that 
remain, the power of rj is 7^2n+2m+i-2n-2fc-i _ ^2m-2k _ ^even power _ sincc the power is 



91 



even and since the number of rj derivatives always exceeds that even power, the derivative 
vanishes,, e.g. 



dr]idr]jdr]k drjidrjj 



V = 4 ^ ^ ^ ^ v'Vm = 4— —^(277,77^ + ri'Smi) = 0. (BIO) 



dr]idr]jdr]kdr]idr]m driidrijdrikdr]i dr]idr]jdr]k 

Note further that the last sum is 

2m+l 



^0 ^ -^0 



(Bll) 



Thus the portion of the integral that "survives" is 



A ' " (2n + 2m + 2)!' ' ^ ' 



and our ten integrals appear as (using V'^rf — l{l — l)rf ^ 



-^1""' ^ i't^ t4 \2m+l/^ TV \2n+1^2n+2m+l 



167r2(2n)!(2m)! 47r(2n + 2m + 2)! 

^^^^^^^ ^ ,(^\ ■ V,)2-+^(^^2 • V,)2'^+2r72"+^-+^ 

(^1 • V^)'"^+'(^2 • V,)'"+'r72"+'"^+^ 

(^1 • V^)'"'+'(^2 • V^)2"+2v2^2n+2m+5 



167r2(2n + 2)!(2m)! 


47r(2n + 2m + 4)!' 




1 


167r2(2n)!(2m + 2)! 


47r(2n + 2m + 4)!' 




1 


167r2(2n + 2)!(2m + 2)! 

-^3 Hill 


47r(2n + 2m + 6)! 



167r2(2n)!(2m + 2)!' 

167r2(2n)!(2m)! "47r(2n + 2m + 2)!^''^''^ ^"^^ ^"^^ ^"^^ 



hrrrn _ 1 . V ,f-^\fj, ■ V ,f-+'r^'-^'-^\ 



167r2(2n)!(2m)! 47r(2n + 2m + 2)! 

— * 

- ^ rV,(^^-V,)^-+^(^^.V,)V"+^"^-^S 



167r2 (2n) ! (2m) ! 47r(2n + 2m + 2) ! 

^8^" \ . V7 ^2m+2/^ . N2n 2n+2m+l 

167r2(2n)!(2m)! " 47r(2n + 2m + 2)! ^'^^ ^''^ ^"^^ ^''^ 

— * 

bm;!!: - \ V if) ■ V l^'^+^ffT . V ^2n+l 2n+2m+3 

167r2(2n)!(2m + 2)!" 47r(2n + 2m + 4)! ^''^'^^ ^''^ ^"^^ 

167r2(2n)!(2m + 2)!"47r(2n + 2m + 4)!^''^'^^ ^"^^ 
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167r2(2n)!(2m)!' ^ ' 

Substitution into Eq.(g§) and Eq.(^ leads to Eq.(|3) and Eq.dOD 
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APPENDIX C: THE EVALUATION OF {ki ■ d)^{K2 ■ d)^ri^+^-^ 
Consider the case with a — 2m + 1, b — 2n + l, 

(«i • df"'+\K2 ■ (ci) 

Let 

K,i = Kiki, K2 = ^2^2 ; k\ = l = k\, (C2) 

and 

^2 ~ OLKi + ; ki ■ = 0; ki- k2 — a ] 
a^ + P^^l; /3 = (l-(«i-«2)')^ (C3) 

so 

Eq.{Cl) = «;^+l«:2n+l 'g' + lV,^2n+l-,^^^ . ^^2m+l+j^^^^ . ^yn+l-j^2m+n+l_ 

j=0 \ J J 

(C4) 

Orient our axes so that 

fj — xi + yj + Zk, (C5) 

with 

Z = f]- {ki X k2), (C6) 

so that, the Z direction is perpendicular to the plane containing ki and k2- Further orient 
axes so that 

f] = xki + yk2 + Zki x k2- (C7) 

Thus 

2^^+! /Or) + l\ 



J=0 



(C8) 
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Consider just the portion 



Let 



Thus 



rj = {p + Z )'2 ; p = X + y = zz* ; z = x + iy, z* = x — iy. (CIO) 



d d dz d dz* d 
dz ~'~ dz* dxdz ^ dx dz*'' ^^^^ 



and 



^Qh+k(^_yn+l-y-kQ2n^+2n+2-h-k(^^^* + (C12) 

Since d^z* = = dz*z, the derivatives become relatively simple 

d':^\zz* + ^2)-+-+^ = {m + n+^y. M^^^. ^ ^2)m+n+i-/.-fc^ 

and 

Q2m+2n+2-h-k ^*h+k ^^^* _j_ ^2'^m+r!.+ i-/i— fe 

2m+2n+2-h-k / cy i i o _ J, _ jL\ 

i=0 V ^ / 

(C14) 

Now 

"^^'^ - + ' ^^^^^ 

Note that the factorial in denominator takes care of cutoff, so we can avoid having to worry 
about upper limit on sum. 

^ (m + n + I - /l - k)\ ^2im+n+l)-h-k-i^^^. ^ ^2y-m-n-f _ .^16) 

(i-m-n-|)! 
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Thus combining factors 

2n+l /o I 1 \ • 2{ni+n+l)+j 2n+l-j (2(m+n+l-/i-fc / „2 ^ w ^* x h+k 

T,{ l Yiy E E E Gf 7 

j=o \ J / y h=o k=o i=o / \ ^ / 

f2m + 1 + j\ f2n+ 1 - j\ f2m + 2n + 2 - h- k\ 1 (n + A;)! 



/i J\ k J\ i J {i — m — n — ^) {h + k — i)\ 



(C17) 



Consider second two lines and let upper limits be oo and let the factorials set the limits. 
That becomes of form 



oo oo oo oo 



EJ2Y.J2fU^h,k,h + k,t). (C18) 

j=0 i=0 h=0 k=0 

Let q = 2n + 1 — j, p = h + k. Then 

/„-^x2n+l oo oo oo p / — '//?\ 9 / r)2 \ V r* \ P 

E,.iClS) = (2„ + 1). I E E £ £ ^ ii, ^ 

11 q\ {2{n + m + 1) - q)\{2{m + m + 1) - p)\ 

X 



g! {2n+l- q)\ k\{q - k)\ {p - k)\(2{m + m + I) - q - p + k)\ 

p\ % 1 

i\[2{m + n + 1) — p — i)\ {p — i)\ {i ~ m — n — ^) ^ ^ 



Use factorials to place upper bound on summation limits. Consider q first; g > 0, q> 
k, q < 2{m + n + l) —p+k;q > 2n + l. Let s = q — k > 0. Upper bound appears ambiguous. 
One can make the replacement 

(2(n + m + l)-g)! _ / d .(n+m+i)-, , _ f d o(n+^+i)-.-. , .p..^ 

and so our s sum is restricted to 

< s < 2(m + n + 1) -p . (C21) 
Note that the right hand side must be positive. So we can perform the g(or s) sum involved 

in (^). It is (using 1 = ) 
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2(m+n+l)-p /_if)\k+s'l 1 

..^2{n+m+l)-k-s ( ± ^ 

a J s\ {2(n + m + I) — p — s)\ 



s=0 



\ a J \ aw J (2(m + n- 



(2(m + n+ 1) -p)\ 

[ a ) r a) (2(m + n+l)-p)!- ^^^^^ 

So our expression (|C19|) reduces to (canceling the last factor) 

r/^yp-yi p! 1 1 1 1 1^^^ ^ (^23) 

,zz* J \zj k\ (p — k)\ i\ (p — i)\ (2{m + n + I) — p — i)\ (i — m — n — ^)\ ^~ 

We can perform the k sum and note that the p sum now has well defined limits since 

2(m + n + l)— p — 1>0 and therefore i < p < 2{m + n + 1) — i. 

Thus 

.,,«3,^(2„.„(-) (-) g g ,.(..J-) 

xL_!^y^'™'YJl)Y£)'_ L ^ ,C24) 

V a J \zz*J \ zj [2[m + n + 1) — p — i)\i\ [i — m — n — ^y. 

Next, noting that w^^l + = (w + and that the p sum produces (2(m+n+i)-2i)! 
see that i is restricted toO<i<m + n + l. 

The p sum is best performed by changing variables to r = p — i. Then < r < 
2{m + n + 1 — i). Thus we obtain 

2n+l / d \ 2m+l rn+n+l 



B,.(C24)^(2„ + l).(^) (A) g 



/3 / \dwj ~^ i\ {2{m + n + 1 — i)\ {i — m — n — ^)\ 



w + iP / a) \ W r] 



2 \ » 



1 + 



z*iw + f) 



\{w — iP/a)J \zz* J \z 
To perform the i sum, notice that 
1 



z(w — i-) 



2{m+n+l-i) / ^^\2(m+n+l) 



w- — ] . (C25) 

a 



{i-m-n- 3/2)! (2(m + n + 1 - i))! 

7r2 



2(™+"+i-*)(m + n + 1 - z)!(m + n + i - z)!(l - m - - 3/2)! 
71^ siniirii — m — n — < _\i-m-n+i 



2' 



22(m+n+i-i) {jn + n + l-i)\ vrh(™+"+i-^)(m + n + 1 - i)! ' 
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(C26) 



Further notice that the 2nd hne of (|C25| ) can be simphfied to 



W -tp /a j yrj ) ^2{m+n+l) 



(C27) 



So 



/-im / rl \ 2m+l ( _\m+n+l i -o 



m+n+1 



■)2i 



E 



7r2 22{'^+"+i) 



Q i! (m + n + 1 — «)! [z{w — ij3/a) + + ijS/ayf''' 



The sum produces 



(C28) 



(m + n + 1)! 



- i[3a) + + 



m+n+l 



(C29) 



Combine this with (|C28| ) and then ( |C17|) gives 



Eg. (CI 7) = /t^'^+i/t^n+i 



\ m+n+1 



v2m+l 



dw' 



7]-2 22(m+n+l) 

(z(w - + z*{w + i(3/a)f - Ar]^{w^ + 



(m + n + 1)! 



(C30) 



Before expanding or simphfying consider case when 



K,2 = = —K, =^ P = 0. 



Then 



Eq.{C30) = (- 



im+n 



^2m+2n+2 

— i (m + n H — )\n 

7ri22(n^+'^+l) 2^ 



-2m— 2n— 3 



(2n+ 1)! X 



d (w2(4a;2 _ ^2))m+n+l 

dwj (m + n + 1)! 



(C31) 



Use 



d'"^^' 2m+2n+i | _ (2(m + n + l))! 

W 10=1 — 



^^2m+l 

Also (m + n + I)! = ^^/r+t+''i)i ' .^r^X.. ■ Thus 
Eq.{C30) = 



(2n+ 1)! 



2m+2n+2 [(2(m + n + 1 ) ) !] ^ 



[(m + n + l)!]222'»+2n+i V r] 
-[(2m + 2n + l)!!]'^(/t' - (/« ■ r))2)™+"+\ 



r^2 \ m+n+1 



(C32) 



(C33) 
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Return to ( |C30|) and consider z{w — ij3 / a) + z* {w + ijd / a) = 2{xw + yP/a). In ( |C30|) this 
factor becomes 



22(m+n+l)(_)m+n+1^^2(^2 _ ^2) ^ a^T]^ - y^) - 2xyW (3 / 0]""^''+^ = 
22(m+„+l)^_^m+n+l[^2^^2 ^ ^2/^2^ _ ^ ^^^^^2p+n+l_ ^(^34) 



Thus 



Ea (C30) - [2(m + n + l)]! a^-^\2n + l)\ 

Eq.[C6U)-K, 22™+2n+2[(^ + ^+l),]2 ^ >< 



UI = 1 



The derivative is of the form 



This radical is 



2xy[3w P'^ f y^ 



with roots 



w = 7 ^ — -. (C38) 



1 - 4 



Since the argument of the root is 



q,2 \^ ^2 ^2 
the roots are complex. Thus 

aw'^ + bw + c = a{w — 7 — — 7 + 

where 



(C35) 



^1 - 4 - 4 ) < 0, (C39) 



a{rf — x'^y a{rf — x'^) 
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Consider the derivative: 

^2m+l 



m+n+1 



J=0 



/^2m + l^ [(m + n + 1)!] 



j=0 



J J (m + n + 1 — j)\(n + j — m)\ 



{w--f- iS)"'+''+^-^{w - 7 + (C41) 
The factorial further restricts the sum to 

min(2m+l,m+n+l) 

E 

j=max(p,n—m) 

Consider the case when (3 — — 5 — ^. Then as w = 1 

mm(2m+l,n+-+l) (^ + ^+1)1 (m + n+l)\ 

Eq.{C41)^ ' ' ' 



j=max{0,n—m) 



j J (m + n + 1) — j)\ (n + j — m)\ 



In order to perform this sum, consider the related product (for s = 1) 



n ./ • ^ (2m + l\ 

= (let / = z — m — n)=2^ . ]^ 

7=0 \ 3 / 1=0 



2«+i /2n + 1^ 



_g£g2n+m 



2(m+n+l) / 
h=0 \ 



But consider the product of the two sums. Let h — j+£so that < h < 2(m + n + l). Let 
I — h — j. Thus, since Z > => h > j and since 2m + 1 > j, we have j < min{2m + l,h). 
Now we also have 2n + 1 > / or j > /i — 2n — 1; thus j > max{0, h — 2n — 1) so 

2(^g+i) /2(m + n + l)\ mM^+i,/.) /2^+iW2^+i\ ^^^^^ 

^=0 V ^ / h=0 j=max{0,h-2n-l) V ^ J \ j J 

Consider the term h = m + n + 1 in the sum. Then 

(2(m + n + l)! _ "'M2mW«+i /2^ + A / 2n + 1 \ 

[(m + n + l)!]2 ^ I j JL+^_^J' ^ ^ 

' / J j=max{0,m—n) \ / \ ' ■J / 
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where we have used = f •) = f )• Thus 

\n—3 J \n+m+l-jj \n+j-mj 

Note that 

Fnir'^^^\ - ,,2(n.+n+i) {2{m + n + l)\f 1/ x^y+n+i _ 
t.q\C6b)y=,--K 22(™+«+i)[(n + m + l)!p;^l^~^y' 

= -([2m + 2n + - . ^)2)™+n+i^ ^Q4g) 

which checks again with Eq.(C33). Thus, with this we write again 
Fn (C-^^) - (2(m + n + l))! / _ 

mm(2m+l,n+m+l) + / 2n + 1 \ 



E 

j=max(0,m— n) 



j I + j — m. 



(w - 7 - - 7 + z5)"+^-'"U=i. (C47) 

To simphfy, let j = k — n. Also, note 



- 7 - i6r+''+^-' {w - 7 + 
1 

^(r^^-x^) -/?(x?/-ir/z)]'=-'". (C4J 



I w=l 
2n+l 



Thus Eg. ( |C47|) becomes 



[(2m + 2n + l)!!]^ I ^ ~ ^ ' mm(2m+n+l,2n+m+l) 



[(m + n+l))!]2 



X 



(2(m + n+l))! r/4"+2 

^ ^1 ^ ^] (^'« - /^(^^ + ^vZ)r'-''''-'-\(r'a - f3{xy - ir^z)f-'"^ . (C49) 

where 

^2 = ry2 _ (^5Q) 



If m > n then sum becomes {I := k — m 

. + 1\ / 2m + 1 
/ j U + m + n 



i=0 
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Because of the symmetry we have 
y-_y- /^2n + lW 2m + 1 ^ 



{a'a - {xy + i7]Z)l3f ''^^-\(T'^a - {xy - ir]Z)PY + c.c. 



(C52) 



The bracket [ ] in ( C52 ) above is of the form 



] = {u- ivf''^^~\u + ivY + c.c. = (let u + iv = re"^) 



^2n+l/gi(2(«-n)-l)0 



+ c.c.) = 272"+^cos(2(n - /) + 1)0, 



(C53) 



in which 



^1 pTlz 



tan 



a'^a — xy(3 



(C54) 



We further need 



2k 



cos{{2{n-l) + l)(j)) = Y. ("^'^ ^"^ ^ 'l(-)^(szr20)2fc(cos(/>^2("-'-^)+^ 



fc=0 
n—l 

E 

k=0 



'2{n-l) + iy 
2k 



cos 



■k+j)+l 



(C55) 



j=0 



where 



cose 



a'^a — xyP 



(C56) 



Recall the original variables 



2 2 I ^72 2 I 2 2 I 2 I ^2 

T] = p + Z = cr + X = X + y + Z , 



Thus 



a = ki- (3 = \Jl- {ki- K2)^, 

X = ki-ff; /3y = ki± ■ rj = {k2 - aki) ■ ff. 

d^a — xy(3 = a{'rf' — (ff - kiY) — fpKi{k2 — aki) ■ ff = 
= rj^i^i ■ k2-ki- 7]k2 ■ fj), 

^22-2^2 ^ ^4j^^ _ . ^^^|2 _ . p.^^, _ . _ 2p.^ . p-^p. . f]k2fi)), 



(C57) 



(C58) 
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and 



Thus 



. Ki ■ K2- Ki ■ r]K2 ■ Tj 
COS(p = , ^ , , 



and 



„2n+l / ™2\ m— n 

^4n+2\ ^2 J 1 2 



= (1 - {k, ■ vrni - ■ vrriKi - {n,r^ff'\4 - {k^ ■ vfY'\ 

So finally factoring out cos4> we obtain 

{ki ■ a)2™+l(K2 • a)2"+lr;2m+2n+l ^ 

= 2^7— TTTT- ^ + m + 1 !p 

(2(m + n + 1))! ^ 

{ki ■k2-ki- f}k2 ■ fj){Kl - {ki ■ fj)f){iil - k2 ■ f})fT 



l=Q k=Oj=0 



I Jyl + m — nj\ 2k J \j 



Check once again the limit ki — — ^2 where Eq.(63) becomes 

where the triple sum J2 collapses to a single sum: 

^(cos0 = 1) =^ j = A; = 0, 

so that 

^ _ ^ f'2n + l\f 2m + 1 \ 1 (2(m + n + 1))! 



^ \ I J\l + m-nJ 2 (n + m + l)!]2 
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APPENDIX D: DERIVATION OF THE INVARIANT MASS M USING THE 

NEW DIRAC BRACKETS 



Kerner has shown [^] that it is possible to develope a single-time Hamiltonian formula- 
tion of Wheeler-Feynman dynamics. His idea, basically is to replace the infinitude of "field" 
coordinates by an infinity of mechnical ones and then with the high order of the equations of 
motion replaced by higher powers of the momentum in the interaction. In Ref . , Crater 
and Yang give a modification of his approach to obtain a Hamiltonian expression for both 
scalar and vector interactions through order The approach taken in this Appendix 



is similar to that given in |Q with two important distinctions: 1) terms of all order in 
1/c^ are included and 2) the effects of the new Dirac brackets are included. The net result 
is an expression for M that agrees exactly with the results obtained in Eg. (|6.14| ). This 



result would not be obtained without the use of these brackets in working our Hamilton's 
equations. 

We begin with the Lagrangian expression for the invariant mass [see Eq. ( |4.10| ) with 
A(r) = 0; hi is defined in Eq.(U)] 

Erel = h{r]i, r]2, r]) = + + — — + QiQ2hi{r]i,r]2, i]) : = 

: = ho{f]i,f]2,f]) + QiQ2hi{f]i,f]2,f]). (Dl) 

In order to find the Hamiltonian H{k,i, k,2] v) from hi{ffi,ff2, ff) we must demand that Hamil- 
ton's equation be satisfied . We use the Dirac bracket since we have used the constraint as 
a strong condition on the dynamical variables. Thus we begin with 

fi, = {ff„Hr = {ffi,H}- 

-[j (f(y{r]i,-^QiA^si{.T,a-f]i,Ki)} ■ {-YlQj^^sji^,^ - Vj,f^j), H} - 

* j 

-{Vi^-J^Qj^^Sji-r,^ - f]j,Kj)} ■ {-J2QiA±siir, ^ - ffi, Ki), H}], (D2) 

j i 

in which 
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H^Ho + QiQ2Hi{Ki, K2]ff), (D3) 

and 



= yjml + kI + ^mi + «i. (D4) 
Substituting this Hamiltonian into the above bracket and using Grassmann truncation yields 



-<5in±si(^, - 571, /^i)} ■ {-Q2A±s2{r, a - 172, ^2), V"^2 + '^i} 



?72 = y + <5l<52- 



But 



^Jml + /?2 ^''^^9k2 
-[ dV{^2, -Q2^±52(t, - ^2, ^2)} • {-Qin±5i(T, a - rfi, Ki), \/ml + kI} - 
-{m, -Q2fi±S2{r, a - 572, K2)} ■ {-<5i^±si(t, a - ffi, ki), ^Jml + fil} = 

^2 ^ Q.q f^^ \ ^ i ^ 1 /C,^ (D6) 

^mi + «i ' '^^2 5«2 y^Jml + ni dm) 



y X Bs)iT, a) = — /C12 = -^^12 := -^^12- 



Hence 



^1 = / . + ^1^2^?^-^ > ■ I dME±s X Bs){T,a) 1 , 



2 
^2 



Substituting this into we find that 
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ho = Ho + — h 

47r?7 



'7712 + 1^2 

^ (m| + ^^^ ■ 4^{Q,Q2H, - , ■ I dME±s X Bs){r, a) 



= H -QiQ2hi = HQ + QiQ2H^-QiQ2hi. (D9) 
Thus letting 

i^i = T^ + ^i, (DIO) 

4:117] 

and using [see Eqs.( |6.4| ) and ( |6.5| ) for the expressions of hi and /ii] 

/^i(rfi,r/2,r/) = /Ii( , 1' ■ ^ jf)+0{QiQ,), (Dll) 



where 



I o?V(^x5 X 55)(r,a) := QiQ2hM,^2,v), (D12) 
then we obtain the following differential equation for Hi 

~ dHi ^ + /?2) 57/^ ^ (^2 ^ g2) ^ 
-ni ■ Ki 2 ■ 2 ~ 

OKi mf 0K2 m2 

, . /?1 /t2 

^mf + /tf -^m^ + K^ 
—2 '^1 ■ ^( / 2 ^2 ■ / 2 ^2 ' / 2 ^2 '"^^ ~ 



"^2 + '^2 V"^! '^i V "^2 + '^2 

—2 ''2-i^(^===-/ii(^=^,^===,r?)+0(QiQ2). (D13) 

The 0{QiQ2) term gives a vanishing contribution as both sides are multiplied by QiQ2- 
Using this and V^|r/|' = /(/ — l)|f/|'~^ we obtain, in addition to the expression given for hi 
in Eq.fCT), 
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— * — * 



QiQ 



oo oo 

2 



Ki K2 

X 



2n+2m+l 



m=on=o yJm{-\-K\ -\- 

V ymj+Kj ym^+K^ ymj+Kj -^"^2+'^2 ' 

(2n + 2m + 2)! 

/ . Ki^mf+K^ ^ \2m+2/ K2 vy '|2n+l 1 / M \2m+l/ K2 vy ^2n+2^ „2n+2m+3 



(2n + 2m + 4)! 



(D14) 



The differential equation for Hi is of the form 



Hi - OHi = hi{ , , , , 77) - d • hi{ , , , rf), (D15) 

where we define the linear operator O by 

= 2 '^l'^"^ 2 '^2--^, (iJlDj 

mf OKi 1712 OK2 

— * 

and the linear operator O by 

In order to solve for Hi we need first to work out the right hand side of its differential 
equation. First note 

Jmi + Ki Jml + nl Jml + Ki Jml + Hl + y^i + ^2 

while 

- [(2n + 1)^4= + (2m + 1)^^=^] x 
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So Eq.(|6.5|) imply 

6 . Q,Q,u^^, ,-, = *^ f f ([(2„ + i)(^ii=^ . V,,/ + 

Wmf + dm2 + K2 m=On=0 dmi + K\ 

+ (2n + 2m + 4)( ■ V,)( ■ V,) + (2m + 1)(^L= ■ V,) 



X 



? + ./m2 + K2 (2n + 2m + 2)! 



mf + nf y m2 + ^2 

Ki VT N2m+1 /^^^2^= Y7 ^Sra+l 

V™i+'^i V'^2+'*2 _^2n+2m+3 



(D20) 



(2n + 2m + 4)! 

In analogy to our decomposition of the field energy integral we decompose this into single 
and double sum pieces giving 

0-QiQ2hi = — — > [^^=^= ■ X 

47r (2m)! ^^2 + ^2 ^^2 + ^2 

, ■ V,)^™ + ( , ■ V,)^™) )r/2-i - 

-^2^lV o , r?2 /„.2 , r?2 



(2m + 2)!V / 2 + ^2 " y^^2 ^ ^2 " 



mf + /?i Y ""^2 + ^2 

00 00 ? ? 

Ki K2 



+ E E ((2n + 2m + 4)[ 



m=on=o ^mf + /5f ym^ + ^2 

V^^+^i A/"^i+«i 2n+2m+l _ 

(2n + 2m + 2)! ' 

r Ki _ y N 2771+2/ K2 . y \2n+2 
^mf+Kj ^m|+Ki 2n+2m+3l , 

(2n + 2m + 4)! ' ^ 

/ Ki . y \2m+2/ K2 . y \2n+2 

+ (2n + 2m + 6)[ , ■ , v/^_! " ' 2^+2^+3 

/ml + k\ ^ml + kI (2n + 2m + 4)! 

y \ 2771+3/ n-2 . y ^2n+3 



^m\+i(\ ^rri^+K^ ^2n+2777+5j ^ j (D21) 



(2n + 2m + 6)! 
Combining like terms we obtain 



V V '^l + ^1 V + ^ V + '^1 V + ^^2 
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^2n+2m+l 



oo oo ? ? 

+ EE(-(2n + 2m + 3)[ , ' , • y ' , X 

m=0 n=0 Y "^1 + /^l ynT'2 + '^ 

(2n + 2m + 2)! 

/ Ki . ^ A2m+2/ K2 . ^ \2n+2 
^mj+Kf ^/rn^+i^l 2n+2m+3l _ 

(2n + 2m + 4)! ' ^ 

-(2n + 2m + 5)[ , ■ , ' 2n+2.+3 

\/mf + K? ^mi + «i (2n + 2m + 4)! 

Ki . ^ A 2m+3 / K2 . ^ ^2«+3 



^ ^' ^2n+2m+5j-^-j (D22) 



(2n + 2m + 6)! 

Based on the above expression we assume a particular solution of the form 
Q1Q2 v^r '^1 '^2 

^" m=on=o y mf + ftf \Jm2 + H2 

( Ki Yj \2m+l( K2 \7 ^2n+l„2n+2m+l 



(2n + 2m + 2)! 

m+2t K2 

(2n + 2m + 4)! 



( Ki Y7 'i2m+2/' K2 VV ^2^^+2 2n+2TO+3 



/■ Ki Y7 N2m+2/ K2 VV \2n+2 2n+2m+3 



Using 



and 



Y^f + K? y^mi + (2n + 2m + 4)! 

(^4 V,,)2"'+-^(^% yj2n+3^^2n+2,„+5 

-cimn^^^S ]_ P23) 

(2n + 2m + 6)! ^ ^ ^ 



0^i= . = 2 . ^i=, (D24) 



'^1 ^4 \2m+l/ '^2 V7^2n+l 



y mf + Y ^2 + 
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= (2n + 2m + 2)( • V^)2-+i(^?= . V^)2"+\ (D25) 

the left hand side of the differential equation for Hi{ki, K2, ff) becomes 
(1 - 0)Q,Q2H,iK,, fi2,f]) = X 

OTT 

i-(^ii=.v,)(^^.v.,J = 

' T7 ^ /___9 , -^9 /___9 , -^9 " / 



'2 
^2 



— — 2^ 2^[-a„„(2n + 2m + 3) • , x 

VT \2m+l/ K2 V7 ^2n+l„2r^+2m+l 



+ 



+6^„(2n + 2m + 3) 



(2n + 2m + 2)! 

^^f:^, ^^^^ ' "'^ 



Ki . j2m+2^ K2 . j2n+2^2n+2m+3 



^2 



(2n + 2m + 4)! 



-c„„(2n + 2m + 5) = • , x 

( Ki vy \2m+2/ K2 VV N2n+2 2n+2m+3 



(2n + 2m + 4)! 



■ry \2m+3/ K2 vy \2n+3 2n+2m+5 

^ V) \ / 2 , r;2 ■ ^'7^ '/ 



Comparing the two sides of the equation leads to 



and thus 

QiQ2-H"i(ki, K2, 77) = - X 



Ki _ yim+l^ K.2 . j2n+1^2n+2m+l 



Q1Q2 ^ /Vj ^V^^^ "'^^ ^ y/^^ 



47r ^0^0 ^m? + «? ^mi + (2n + 2m + 2)! 

Ki vy \2m+2/' K2 vy \2n+2 2n+2m+3 



(2n + 2m + 4)! 



2 

+ 
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( Ki \2m+2( K2 VT \2n+2 2n+2m+3 



ml + nj \lml + kI (2n + 2m + 4)! 

K2 Y7 \2n+3„2n+2m+5 

(D28) 



N2m+3/ K2 VT \2n+3 2n+2m+5 



(2n + 2m + 6)! " 
which agrees exactly with the Darwin portion of M obtained earher Eq. (|6.14|) . 



Ill 



APPENDIX E: SCHILD-LIKE SOLUTION FOR THE TWO-BODY PROBLEM IN 

THE CASE OF EQUAL MASSES 



Here we present the semiclassical Hamilton equations for the equal mass case restricted 
to circular orbits (we suppress the tilde notation). First note that the Hamiltonian Eq.( p.37D 



is of the form 



(El) 



H = 2vW^ + ^ + ^/(^^ (/? ■ vf). (E2) 

4:717] Stti] 



Thus, Hamilton's equations are = df /Ok^, f^2 = df /d{K ■ fj) 
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-li = ^ = —. — jri - - — jrif{K ,{K-r]) ) + /,2 [k ,{K-r]) ) k{k ■ r]), 

From this we can see that circular orbits defined hj ff ■ ff = are implied hj k, ■ ff = 0. This 
furthermore implies that —k, ■ k, = 0. Thus not only is ff = const, but also = const. 
Imposing these conditions on the above Hamilton equations we can simplify our equations 

above to 



^ ^ ^ _ ^1^2 -. Q1Q2 ^jrp2 QN _ 



OK ^J K? + 4:7ir] 
with B and A constants. Combine the two equations and we find 

f^=J = -BV, (E5) 

so that 

f] = -ABf]:=-n'^f], (E6) 
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with 



Q1Q2 ri ^2 2/?V^' + m2 



m 



The frequency is defined by the initial data rj and \k\ and is real for Q1Q2 < for < < 
in which is the value at which = 0. Let us remark that at the semiclassical 



level Qi are Grassmann variables: therefore ~ \fQ\Q2 is to be interpreted as an even 
algebraic object satisfying fi^ = 0, = Q2^^ = 0. 

We also find 



H = 2V/t2 + m2[l - -n'^ff]. (E8) 



Our solution to Eg. (|E6|) is 



and thus 



r/(r) = a cos Qt + f3 sin Qt = 

= 5(1 - ifiV^) + /3(fir - ^n'r'), (E9) 



^(r) = -aQ\ + - -l^V^). (ElO) 



But rf ■ ff = 0. This implies that 

- + a ■ - 2nV2) = 0, 

so that 

32 ->2 



(3' = a\ a-/? = 0, (Ell) 



and therefore 



a" = ff. (E12) 

Otherwise the vectors a and (3 are arbitrary. 
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